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Abstract 

The scaling theory of the transitions between plateaus of the Hall conduc- 
tivity in the integer Quantum Hall effect is reviewed. In the model of two- 
dimensional noninteracting electrons in strong magnetic fields the transitions 
are disorder-induced localization-delocalization transitions. While experimen- 
tal and analytical approaches are surveyed, the main emphasis is on numerical 
studies, which successfully describe the experiments. The theoretical models 
for disordered systems are described in detail. An overview of the finite- 
size scaling theory and its relation to Anderson localization is given. The 
field-theoretical approach to the localization problem is outlined. Numerical 
methods for the calculation of scaling quantities, in particular the localization 
length, are detailed. The properties of local observables at the localization- 
delocalization transition are discussed in terms of multifractal measures. Fi- 
nally, the results of extensive numerical investigations are compared with 
experimental findings. 
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I. INTRODUCTION 

As announced in a paper by von Klitzing, Dorda, and Pepper (1980), Klaus von Kl- 
itzing discovered that under certain experimental conditions the Hall conductivity a^y oi 
a quasi-two-dimensional electron gas is quantized to a high precision in integer multiples 
of e^/h while at the same time the longitudinal conductivity a^x vanishes. A remarkable 
aspect of this integer quantum Hall effect (QHE) is that the quantization exists over a finite 
range of physical parameters, like the magnetic field B or the carrier concentration n. This 
observation cannot be explained by the semiclassical Drude theory of conductivity. If the 
longitudinal conductivity cr^.^. vanishes the Drude result for the Hall conductivity is given 
by (Txy = ve^/h, where v = nh/eB = n2Til1 is the filling factor. Only if v is an integer, 
i.e., if the Fermi level lies exactly between two Landau levels, is the Hall conductivity quan- 
tized. Laughlin (1981), Aoki and Ando (1981), and Halperin (1982) showed that the Hall 
conductivity takes on its quantized value ie^/h over a finite range in filling factor around 
z/ = i if in that range the states at the Fermi energy are localized. By the same argument, 
when the Hall conductivity is quantized there are states helow the Fermi energy that are 
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not localized. The width of the quantized plateaus then depends on the ratio of localized 
vs. extended states. 

The concept of localized states in disordered systems was developed by Anderson (1958). 
He showed that if a quantum-mechanical system is sufficiently disordered, states have a finite 
probability to return to a given site in the long-time limit. This absence of diffusion implies 
that these states are localized in a finite region of space. The transmission probability decays 
exponentially on a length scale, which is called the localization length. For locahzed states 
the static conductivity vanishes at zero temperature. On the other hand, if the disorder 
is weak enough extended states might exist that do not decay exponentially and fill the 
whole system. Their contribution to the conductivity is finite even at zero temperature. 
The energy that separates extended from locahzed states is called the mobility edge. At 
the mobility edge the character of eigenstates is different from both extended and localized 
states and we will call these states critical states. 

The observation that an explanation of the QHE involves both extended and localized 
states was the more unanticipated as the scaling theory of localization (Wegner, 1976; Abra- 
hams et al., 1979; Wegner, 1979) predicted the absence of extended states in two-dimensional 
systems. It is the presence of a strong magnetic field that leads to the emergence of nonlo- 
calized states in two dimensions (Ono, 1982b). Chalker (1987) showed that these states exist 
only at a single energy in the limit that scattering between Landau levels can be neglected. 
At zero temperature the Hall conductivity is thus expected to exhibit sharp steps whenever 
the Fermi energy passes the critical energy. The longitudinal conductivity vanishes for all 
energies except at the critical energy. The quantum Hall plateau transition is thus a spe- 
cial case of a metal-insulator transition. The insulating phases correspond to the plateau 
regions, in which a^x vanishes and a,j.y is quantized. These insulating phases are separated 
by mobility edges with a finite Uxx- Truly metallic phases are absent in the quantum Hall 
system. The states at the mobility edge show critical fluctuations (Wegner, 1980; Aoki, 
1983; Aoki, 1986). 

The description of this localization-delocalization transition in terms of a zero- 
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temperature phase transition is the topic of this review. We will focus on the theory of 
noninteracting electrons. Thus we do not consider the fractional QHE (Tsui et ai, 1982), 
where the quantization of the Hall conductivity in rational fractions of e^/ h is due to electron- 
electron interactions. One notes, however, that experimentally the fractional quantum Hall 
transitions show a remarkable similarity to the integer transitions (see Sec. y). In the ab- 
sence of an analytical theory that is capable of providing quantitative results we will mostly 
rely on numerical simulations. 

The review is organized as follows: In Sec. |T| we shortly survey the experimental in- 
vestigations of the quantum Hall transitions. Section |lll| discusses the different models of 
disorder used in theoretical calculations. Section |^ provides an overview of the finite-size 
scaling theory. Section outlines the description of the Anderson transition in terms of 
a field theory and its extension to the quantum Hall system. Sections and |VII| explain 



which physical quantities are used in numerical finite-size scaling studies and how they are 
computed. In Sec. |V1II| the results of numerical calculations are presented and analyzed in 



terms of the finite-size scaling theory. In Sec. ^ we characterize the critical states using a 
multifractal analysis. In Sec. ^ we compare the experimental and theoretical results pre- 
sented in the previous sections and discuss their implications for an understanding of the 
quantum Hall transition. Except where noted otherwise we try to provide enough detail to 
allow the reader to use this review as a starting point for own investigations of the subject. 

A comprehensive introduction to the QHE is provided in the book edited by Prange and 
Girvin (1987). For reviews about disordered electronic systems and localization see Thouless 
(1974), Lee and Ramakrishnan (1985), and Kramer and MacKinnon (1993). 

II. EXPERIMENTAL RESULTS 

In this section we want to review experimental investigations of the transition between 
plateaus in the Hall conductivity. We will focus on those results that are uniquely related 
to these phase transitions. In later sections we will show how a simple scaling theory for 
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the integer quantum Hall transitions accounts for most of the observed features. We will 
therefore postpone the discussion of the experimental results until we are familiar with the 
necessary theoretical framework (Sec. 

The sections [11 A] to [11 C| are ordered roughly chronologically. We will start with exper- 
iments on the temperature dependence of the transition. Although these experiments give 
the most indirect information about the scaling behavior they were the first to be performed 
and instrumental in developing the idea that the plateau transitions are continuous phase 
transitions. Much more direct contact with the concept of finite-size scaling is made in 
experiments with varying sample sizes. Finally, the dynamic aspects of the transition are 
elucidated in high-frequency experiments. 

A. Temperature-Dependent Scaling 

Soon after the discovery of the QHE (von Klitzing et ai, 1980) it became apparent 
that the quantized plateaus can become extremely broad and the transition between them 
extremely sharp as the temperature was lowered (Kawaji and Wakabayashi, 1981; Paalanen 
et ai, 1982). At the same time the longitudinal conductivity exhibits a series of sharp spikes 
at the position of the plateau transitions. Since the Fermi energy is in a region of localized 
states when the Hall conductivity is quantized and the longitudinal conductivity vanishes, 
it can be concluded from the wide plateaus that most of the electronic states are localized at 
low temperatures. Paalanen, Tsui, and Gossard (1982) estimated that in an AlGaAs/GaAs 
heterostructure at 50 mK 95% of the states in each Landau level are localized. 

The temperature dependences of the longitudinal conductivity a^^ and the Hall conduc- 
tivity of a Si-MOSFET are shown in Fig. |l[ Even in the absence of localization the 
temperature dependence of the Fermi-Dirac distribution leads to a temperature dependence 
of the conductivity, 

a,AT) = JdE (-^) ^,AT = 0). (2.1) 
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If (Jxx{T = 0) is finite only in an interval of width A near the center of each Landau band, 
Eq. ( |2.1| ) predicts an increase of a^x in the band center with decreasing temperature as long 
as fc^T ^ A. This increase is seen in the peaks of Fig. |1|. 

When plotting a^x versus axy as a function of temperature and filling factor striking 
features appear. The data for a Si-MOSFET for temperatures between 0.35 K and about 1 K 
in Fig. 1^ lie on a single curve indicating that axx and axy are not independent but dependent 
on a single parameter. At higher temperatures the effect of the Fermi distribution is still 
visible in the increase of axx- By contrast, the conductivities seem to be independent and 
follow a two-parameter scaling flow in the temperature driven flow diagram for an InGaAs- 
InP heterostructure in Fig. |^. In both figures a flow towards {a^^jaxy) = (0, integer)e^//i, 
corresponding to the widening plateaus, is obvious. Kravchenko and Pudalov (1989) and 
Kravchenko et al. (1990) studied Si-MOSFET and AlGaAs/GaAs heterostructures and 
observed both types of behavior depending on the filling factor. A single curve was also 
obtained by McEuen et al. (1990) in a study of the non-local transport properties in high- 
mobility GaAs-AlGaAs samples. 

It is experimentally extremely difficult to obtain reliable flow-diagrams close to the tran- 
sitions at half-integer Hall conductivity. As the temperature approaches zero the transition 
regions become very narrow and inhomogeneities in the electron density will let the transi- 
tion happen at slightly different values of the magnetic field in different parts of the sample. 
Since axx and axy are measured on different parts of the sample this will lead to irregular 
behavior near the transition. Wei, Tsui, Paalanen, and Pruisken (1988) avoided this problem 
by studying the critical behavior of each component of the resistivity tensor separately. Fig. ^ 
shows the transport coefficients pxy, Pxx, and the derivative dpxy/dB for an InGaAs-InP het- 
erostructure. Note the similarity between p^x and dpxy/dB. However, while the maximum 
values of the former decrease, the maximum values of the latter increase with decreasing 
temperature. In Fig. |^ the temperature dependence of the maximum {dpxy/dB)^'^^ and the 
inverse width (AB)^^ of the pxx peaks is plotted for several Landau levels. Over the range 
of temperatures from 0.1 K to 4.2 K they show power-law behavior {dpxy/ dB)^^^ oc T"*^ and 
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AB oc T'^ with k, = 0.42 ± 0.04. Further measurements of higher-order derivatives showed 
that the extrema of (Pp^x/dB"^ and d^p^y/dB^ diverge hke T"^** and T~^'^, respectively (Wei 
et ai, 1990). In experiments on low-mobihty AlGaAs/GaAs samples the same exponent was 
observed but the characteristic temperature below which this scaling was observed was much 
lower (200 mK) than in the InGaAs/InP samples (4.2 K) (Wei et al, 1992). Huckestein et 
ai (1991) studied the temperature dependence of the plateau width in AlGaAs/GaAs het- 
erostructures and found agreement with k = 0.42 for temperatures below 4.2 K. For Landau 
levels where the spin splitting was not resolved (dp^y/dB)^^^ and (AB)^^ diverge like T~'^/^ 
(Wei et ai, 1990; Hwang et ai, 1993). The same exponent n was measured in the fractional 
QHE regime of high-mobility AlGaAs/GaAs heterostructures for the scaling between filling 
factors 2/5 and 1/3 (Engel et al, 1990). 

The universality of the exponent k, was questioned by Wakabayashi et al. (1989; 1990) 
and by Koch et al. (1991a). Wakabayashi et al. studied two of the valley- and spin-split 
subbands of Landau levels n = and n = 1 in a Si-MOSFET. For the rather limited 
temperature range between 0.35 K and 1 K they observe k, = 0.29 ±0.10 and n = 0.16 ±0.02 
for subbands (0 | — ) and (1 t ~), respectively. Below 0.2 K the temperature dependence 
saturates, presumably due to Joule heating (Wakabayashi et al, 1990). Koch et al. found 
a dependence of n on the mobility of AlGaAs/GaAs heterostructures deliberately doped to 
decrease the mobility. They observed power-law scaling with temperature over the range 40 
mK to 1.1 K. For spin-split Landau levels the measured values for k increased from 0.28 to 
0.81 with decreasing mobility. If the spin splitting is not resolved, k is considerably smaller 
but not by a universal factor of two as found by Wei et al. (Koch et ai, 1991a). 

In the context of the finite-size scaling theory (see Sec. pV]) the power law behavior of 
the transport coefficients reflects the single-parameter scaling of the conductivity tensor 
(Pruisken, 1988) 



where B* is the critical magnetic field of the plateau transition, L^s is the effective system 




(2.2) 
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size, and v is the critical exponent of the locahzation length [cf. Eq. ( [4.3|) ]. The strongest 
support for this assertion comes from the peculiar temperature dependence of the magnetic 
field derivatives of the conductivity tensor. At T = 0, Leg is equal to the system size. At finite 
temperatures the effective system size is given by the phase coherence length L$ (Thouless, 
1977). If L<j, diverges as T"^/^ as the temperature approaches zero then l}^^ oc T~'^ with 
K = p/2v. With Eq. (|2.2| ) the n-th derivative of the conductivity tensor at the critical point 



IS 



unless it vanishes due to symmetry. 



oc oc T-"", (2.3) 



B. Size-Dependent Scaling 

Koch et al. (1991b) measured the exponent v directly by studying samples of the same 
shape but different sizes. For sufficiently small samples (dp^y/dB)^^^ and AB saturate at 
low temperatures (see Fig. |^). The saturation temperature decreases with increasing system 
size. This is interpreted as the temperature where the phase coherence length L$ becomes 
comparable to the system size and the temperature-dependent scaling at higher temperature 
crosses over to size-dependent scaling. The saturation value of AB at low temperatures is 
then given by the condition L/^{AB) ^ 1, i.e., AB oc L^^/^ . By fitting their data to this 
relation Koch et al. obtain v = 2.3 ± 0.1 for the three lowest Landau levels. From the 
measured values of k they deduce that p is not universal and varies between 2.7 and 3.4. 
For a Landau level in which the spin splitting was not resolved u = 6.5 ± 0.6 was measured 
while the value of p was comparable to its value in other Landau levels. 



C. Frequency-Dependent Scaling 



The dynamical conductivity (Jxx{f) of a low-mobility AlGaAs/GaAs heterostructure was 
investigated by Engel et al. (1993) by measuring the attenuation of a coplanar transmis- 
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sion line on the sample surface. The losses in the transmission line are due to the finite 
conductivity of the two-dimensional electron gas below the surface. In contrast to previous 
measurements using crossed rectangular waveguides (Kuchar et ai, 1986; Huckestein et al, 
1991) that were limited to a single frequency / = 35 GHz Engel et al. were able to sweep 
the frequency from 0.2 to 14 GHz. As can be seen from Fig. |^ the effects of raising the 
temperature and increasing the frequency are quite similar. With decreasing frequency the 
width AB of the peaks in Re{axx) decreases and saturates below a temperature-dependent 
frequency. Above this frequency Ai? scales like 

AB oc f\ (2.4) 

with 7 = 0.41 ± 0.04 for spin-split Landau levels and 7 = 0.21 for non-split Landau levels. 
The saturation temperature corresponds to hf ~ ksT in agreement with the results of 
Huckestein et al. (1991). 

Similar to Eq. ( |2.2| ) a dynamical scaling ansatz for (Txx{L,f) is 

aUL, f) = ^SxxiL/^B), from, (2.5) 
h 

where to{B) oc ^{B)^ oc \B ~ B*\^'^^ is the relaxation time of the system that diverges as the 
critical point is approached. For sufficiently high frequencies fTo{B) dominates the scaling 
behavior of S^x and AB scales like f^/'^^. Equation ( p.4|) thus implies 7 = l/z/2;, and using 
the result u = 2.3 from the size-dependent scaling experiments, z = 1 is obtained (Engel 
et al, 1993). 

III. MODELS OF DISORDER 

Before we can attempt to understand the experiments we have to define our system in 
theoretical terms. We will make several simplifying assumptions: (1) The electron system is 
strictly two-dimensional. (2) The interactions between the electrons can be neglected. (3) 
The spin degrees of freedom can be neglected. (4) Boundary effects, like edge states, are not 
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essential in understanding the critical behavior of the system. The electron gas in real devices 
is only quasi-two-dimensional, in the sense that the motion perpendicular to the plane of 
the electron gas is quantized into electrical subbands. If only the lowest electrical subband 
is occupied the electrons have only two spatial degrees of freedom left and we may treat 
them as two-dimensional, justifying our assumption (1). Assumption (2) lacks an a priori 
justification. The treatment of non-interacting electrons will turn out to be quite successful 
in explaining the experimental static-scaling results. While this might serve as an a posteriori 
justification of our approach a theoretical assessment of the relevance of electron-electron 
interactions is still missing. For an understanding of the observed dynamical scaling it seems 
to be essential to go beyond the non-interacting electron approximation. Assumption (3) 
seems justified since most experiments deal with completely spin-split Landau level and we 
neglect interactions between the electrons. For non-spin-split levels spin-orbit scattering 
might be important and one would have to go beyond the spinless-electron approximation. 
The last assumption (4) treats the phase transition in the spirit of standard scahng theory 
as a bulk critical phenomenon. We will make use of this assumption by choosing boundary 
conditions most suitable to the calculation at hand. 

We will thus describe our system by a single-particle Hamiltonian 

H^Ho + Vir), (3.1) 

where 

Ho^:^{p-eAy (3.2) 

is the kinetic energy of the electron and V{r) is the disorder potential. In the absence of 
disorder, V{r) = 0, the spectrum of H for an infinite system is a set equidistant levels, 

En = {n+ l/2)?iujc, (3.3) 

separated by the cyclotron energy hujc = fieB/m. Each Landau level n is infinitely degen- 
erate with a degeneracy of 
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UB = l/27rZ^ (3.4) 

states per unit area. The length Ic — (^/eS)^/^ is called the magnetic length and is the 
classical cyclotron radius in the lowest Landau level n — Q. 

In order to fully define the model system we will now describe different models of the 
disorder potential typically used in calculations. 



A. Real-Space Models 

The statistical properties of a disorder potential are completely determined by the joint 
probability distribution P\y\ of the disorder potential V{y) for all coordinates r. Starting 
point for most analytical calculations is the white-noise potential distribution 

P[V] = AAexp I dh[V{v)f^ , (3.5) 

where J\f normalizes the distribution. With this distribution the potential at different coor- 
dinates is uncorrelated, 

V{r)V{r')=V,'d{r-r'), (3.6) 

and has a Gaussian distribution at each point in space. The overbar denotes the average 
with respect to the distribution P[V] of the disorder potential. This potential distribution is 
well suited for analytic calculations since it allows one to perform the average over disorder, 
due to the Gaussian dependence on the potential. For the white-noise distribution Wegner 
(1983) showed that the density of states in the lowest Landau level is given by 



where 




(3.7a) 



(3.7b) 



12 



If the integrand in the exponent of Eq. ( |3.5| ) contains a more comphcated function of 
V{r), but still depends only on the potential at one point in space, the potential is still 
5-correlated, but higher-order cumulants of the potential at the same site have non-zero 
coefficients. Brezin, Gross, and Itzykson (1984) generalized Wegner's result ( |3.7| ) to this 
more general class of potentials. 

For numerical calculations it is not practical to specify a random value for the potential 
at every point in space. Instead one uses an approximation to the white-noise distribution 
Eq. ( |3.5|) , the (5-scatterers potential, 

N 

V{r)=Y^V,6{r-r,), (3.8) 

where the sites rj of the scatterers are randomly chosen and there is an equal number of 
attractive Vi = —V and repulsive scatterers Vi = +V. This potential has the correlation 
function (p.6|) and approaches the white-noise potential ( p.5|) with Vq = riiV'^ in the limit 
that the density of scattering sites rij becomes infinite. Due to the symmetry of the potential 
the density of states is symmetric with respect to the center of the Landau level. The 
self-consistent Born approximation (SCBA) for the density of states of this potential is a 
semi-circle, 

with the bandwidth T = {AmV^ /2TTliy/^ = 2Vb/\/27r/c (Ando and Uemura, 1974). Since 
this is also the energy scale in Eq. ( p.7|) we will frequently use F as a measure of the width 
of the Landau level. 

In numerical calculations care has to be taken that the concentration q = 27rl'^ni of 
scatterers is large enough in order to avoid a singularity in the density of states due to 
wavefunctions avoiding the scatterers (Ando and Aoki, 1985). To see the origin of this effect 
consider a square of side L containing an integral number of flux quanta Ns = L^/27rZ^. 
From Eq. ( ^.41 ) we see that Ns is also the number of degenerate eigenstates of the system 
in the absence of disorder. We can form a superposition of these Ns states that satisfies Ns 
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(3.9) 



constraints. Thus if the number Ni = riill of (5-scatterers is less than A^^, the wavefunction 
can be made to vanish at the site of each scatterer and hence the energy of this state will 
not be influenced by the disorder potential. This leads to a divergence of the density of 
states sX E = huJc/2 even in the presence of disorder. In terms of the concentration q, this 
divergence will occur for Cj < 1. Ando and Aoki state that for Cj = 40 a good approximation 
to the white- noise limit is achieved (Ando and Aoki, 1985). 

Disorder potentials with a finite correlation length can be generated by replacing the 



5-functions in Eq. (|3.8| ) by functions with a finite range. For Gaussian scatterers, 

Ns 

V{v) = Y.v{v-v,), (3.10) 



1 



v{v) OC r^e^l'^l'/''^', (3.11) 



the correlation function becomes Gaussian, 



1 



VivWiv') oc ^e-l'-'-'l'/^-'. (3.12) 



The concentration necessary to reach the high concentration limit is much lower for Gaussian 
scatterers than for (5-scatterers. The number of impurities an electron can feel increases with 
the range of the scatterers so that for Gaussian scatterers the effective impurity concentration 
is enhanced by a factor 1+cr^//^ relative to (5-scatterers (Ando and Aoki, 1985). For scatterers 
of the same strength V the bandwidth is reduced by a factor of /5 = (1 + a"^ /llY^"^ relative 
to the 5-scatterers, Vfj = V/ (3 (Ando and Uemura, 1974). 

Choosing different numbers of attractive and repulsive scatterers allows for the study of 
potentials that are not particle-hole symmetric on average. A potential with only repulsive 
scatterers was used by Huo et al. (1993) to show that in this case the critical energy is not 
simply related to features in the density of states. 



B. Landau-Space Models 



This section gives an overview over models that specify the disorder distribution in the 
space of the Landau states. While it provides the reader with the basic ideas and results a 
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more detailed description is given in the Appendix 

In a strong magnetic field it is helpful to choose a set of basis states in which the kinetic 
energy part of the Hamiltonian Hq ( |3.2| ) is diagonal. In the Landau gauge 

A = Bxey (3.13) 

and for a strip geometry of width Ly with periodic boundary conditions in the y-direction 
these states are the Landau states 

1 / r — kP\ 

^„fc(r) = {r\nk) = ^=e^'^^Xn , (3.14) 



where 

Xnix) = {2^n\^y^'^ H^{x)e-^"/^ (3.15) 

are the harmonic-oscillator eigenfunctions, and Hn{x) are Hermite polynomials. In terms of 
these states the Hamiltonian including the disorder potential V{r) can be written as 

H = J2Y.\^'^)(^f^\H\nk'){n'k'\, (3.16) 

nk n'k' 

{nk\H\nk') = (n + - ) huJcSn,n'Sk,k' + {nk\V\nk'). (3-17) 



2. 

For a particular realization of the real-space potential, the matrix elements (nk\V\n'k') can 
now be calculated. Instead, one could try to calculate the distribution function of the 
matrix elements given a particular distribution of the real space potential. However, to 
study the universal critical behavior of the Hamiltonian it is not necessary to know the 
whole distribution. The universal properties presumably only depend on a few statistical 
properties of the Hamiltonian, like the correlation function 



f{r,r') = V{r)V{r'). (3.18) 

Under this assumption it is not necessary to construct the entire distribution function of the 
matrix elements. By the same argument only the second moment of the matrix elements is 
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important for the critical behavior. Thus it suffices to construct a random matrix {nk\V\n'k') 
[which we will call a random Landau matrix (RLM)] with a second moment corresponding 
to ( p.l8| ). The correlation function for the matrix elements can be calculated from the 
real-space correlation function /(r,r') 



{niki\V\n2k2){nzk^\V\nik^) 

cfrcfr'{niki\r){r\n2k2) {ri'^k^lr') (r'|n4/i;4) V"(r)y(r'), 



rf'^^'^Kfci(r)^n.fe(r)C3,3(r')^„,fc,(r')/(r,r'). (3.19) 
If the real-space correlation function is translationally invariant, 

f{r,r')=g{\r-r'\), (3.20) 

Eq. ( |3.19D simplifies to 

{niki\V\n2k2){n3k3\V\n4ki) = -^j^Sk^^kiM-ks J dp^^giPx, h - h) J dX 



y c 



[ X + pj2-kill \ fx + pj2-k2{ 

^Xni J Xn2 



c 



( X-pj2-k,ll \ ( X-pj2-k,ll \ 
xxn3(^ 1 )^-\ 1 j' (3.21) 

where g{x^ k) is the Fourier transform of g{x,y) with respect to y, 

g{x,k)= f"^"^' dyg{x,y)e'''y. (3.22) 

J-Ly/2 

The (5-function in Eq. (|3.21|) reflects the translational invariance of the correlation function. 
Equation ( p.21| ) describes the statistical properties of the RLM. For magnetic fields much 



stronger than the disorder potential the coupling between the states within one Landau 
level becomes much stronger than the coupling between states in different Landau levels. 
Neglecting the coupling between different Landau levels leads to the single-band approxi- 
mation in which the RLM becomes block diagonal in the Landau level indices n^. In this 
approximation the physics of each Landau level can be discussed independently. For simplic- 
ity, let us restrict our discussion of the statistical properties of the RLM in the single-band 
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approximation to the lowest Landau level n = and the case of a random potential with 
Gaussian correlations 



The correlation function (|3.21|) then becomes (up to terms vanishing exponentially in the 
limit Ly oo) 

(0fci|\/|0fc2)(0fc3|V-|0fc4) = 6k,-k,M-ksexp (JAz^hl\ JdxJ dx' 



( {X-X'f\ (X-HW (x-k2ll\ (x'-k^ll\ fx'-kj^; 

X exp , Xo ; Xo ; Xo ; ] Xo ' 



^0 



'21x1 iTp'"''^ " ^'^'^'^') ^''P ~ ^'^'^'^) ^'^-'^''^-'^^ ^3.24) 

where jS"^ = {(y'^+lD/l'l characterizes the correlation length of the potential projected onto the 
lowest Landau level relative to the magnetic length. Due to the Kronecker 5 in Eq. (p.24|) 



the matrix elements (0fci|V"|0fc2) are correlated only if they belong to the same diagonal 
ki — k2 = const. More precisely, if (0A;i|l^|0fc2) is split into its real and imaginary parts, 
(0/ci|1/|0fc2)R and {0ki\V\0k2)i, respectively, 

{0ki\V\0k2) = (0A;i|V|0A;2)r + i{0k,\V\0k2h, (3.25) 

then, due to hermiticity, 

{0ki\V\0k2)R{0ks\V\0ki)n = ^{6k,-k,M-k, + Sk.-k.M-kM^, 0, 0, 0; k^, k2, k^, k^), (3.26) 

1 



{0ki\V\0k2)i{0k3\V\0k4)i = -{dk^-k^M-kz - hi-k2M-ki)p{^, 0, 0, 0; fci, k2, ks, k^), (3.27) 

where p(0, 0, 0, 0; fci, A;2, A^a, k^) is the correlation function of Eq. (|3.24]) without the Kronecker 
5, and real and imaginary parts are uncorrelated (Huckestein and Kramer, 1989). The first 
exponential factor in Eq. ( p.24| ) leads to a band structure in the RLM, i.e., as ki — k2 becomes 



large the average magnitude of the matrix element |(0fci|^|0fc2)p becomes exponentially 
small. The second exponential factor governs the strength of the correlations along the 
diagonals. Note that even for a 0, i.e., for 5-correlated potential, the range of the 
correlations does not vanish due to the Gaussian spread of the Landau wavef unctions. 
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We now want to express the correlated matrix elements {niki\V\n2k2) , which satisfy the 
correlation function ( p.21| ), in terms of uncorrelated random variables, that can relatively 



easily be generated on a computer. The correlation function motivates the ansatz 

{niki\V\n2k2) = -—- f dxdx' h{x — x' , ki — k2)uo{x', ki — ^2) 

( X — kill\ ( X — k2ll\ 

where the Uq{x, k) are uncorrelated, complex random variables satisfying 



Mo(a;, k)uo{x', k') = 6{x — x')6k-k'- (3.29) 

For the simple case of Gaussian correlations the intra-Landau-level matrix elements in the 
lowest Landau level are given by 

(0.,|V|0.2)-^exp^ 1 j 



X jdiuApi+ (fci - fc2)/.) exp(-e^). (3.30) 



Note that this is also the form of Eq. ( [Aq ) for any Landau level iii = 712 in the limit 
a /Ic 00 (Huckestein, 1992). 



It is worthwhile to compare the RLM model with the tight-binding model of section |III C. 



The correlation function in ( p.21D describes a random Hamiltonian on a one- dimensional lat- 



tice where the lattice sites are labeled by the integers i = kLy/2'K. The diagonal matrix 
elements {nk\V\nk) correspond to the random site energies e^. The differences to the An- 
derson Hamiltonian are in the off-diagonal matrix elements that correspond to the hopping 
matrix elements Vij. While in the Anderson case the off-diagonal matrix elements are real 
constants^ and couple only nearest neighbors, they are random complex variables in the 
RLM that couple all sites within the range |i — j | ~ Ly/2TTlcf3. As is shown in section |V111 



^In one dimension the hopping matrix elements can be made real by a global gauge transformation 
since closed paths can enclose no flux. 
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for finite Ly, the RLM has only localized states. However, the largest localization length is 
of the order of Ly and in the two-dimensional limit, Ly oo, both the range of the hopping 
matrix elements and the maximum localization length diverge. 

C. Lattice Models 

The standard model of localization is Anderson's tight-binding Hamiltonian (Anderson, 
1958) 

H = J2eS){t\ + Y.V.^^){3l (3.31) 

« {id} 

where i labels the sites on a lattice and are nearest neighbors. The site energies are 

independent random variables and in the absence of a magnetic field the hopping matrix 
elements are constant, Vij = V, and define the energy scale. With this Hamiltonian the 
Schrodinger equation takes the form 

Sitti + Vijaj = Eai, (3.32) 

where a, is the amplitudes of the wave function at site i. For = the Hamiltonian 
describes an electron moving in a perfect crystal. 

This model can be understood as the simplest form of a more general tight binding model 
where the Vij are a function of i and j. Such a model can be derived from the Schrodinger 
equation for a free particle in a periodic potential by expanding the wavefunction in terms 
of a complete set of orthonormal Wannier functions localized at each atomic site (Thouless, 
1974) 

^(r) = EE4"V("nr-R.), (3.33) 

i n 

where R, is the position of atom i and n labels the atomic orbitals. Taking into account 
only a single orbital per site leads to an equation of the form ( |3.32| ). 
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In two dimensions on a square lattice, Eq. ( p.32|) has eigenvalues E in the band between 



—AV and +AV. In the presence of disorder all states in the system are localized (Abrahams 
et ai, 1979; Wegner, 1979; MacKinnon and Kramer, 1981). 



When a magnetic field is applied the Hamiltonian ( |3.31| ) has to be modified. To lowest 
order in the magnetic field the infiuence of the magnetic field on the Wannier functions can 
be neglected and the effect of the magnetic field is to change the phase of the wavefunction 
between two sites. This is achieved by replacing the constant hopping matrix elements by 

Vi, = Vexp [-J rfrA(r)) , (3.34) 

a procedure known as Peierls substitution (Peierls, 1933; Luttinger, 1951). In the absence of 
disorder this system shows a very rich band structure, the self-similar "Hofstadter butterfly" 
(AzbeF, 1964; Hofstadter, 1976). In particular, if the number of flux quanta per unit cell 
of the lattice, a^/27r/^, where a is the lattice constant, is a rational number p/g, the tight 
binding band splits into q magnetic subbands. 

In the presence of disorder this model can quite accurately describe the lowest Landau 
levels of a continuum model. This is the case in an intermediate range of the disorder where 
the disorder broadening of the magnetic subbands is large compared to their intrinsic width 
but small compared to the separation of the magnetic subbands. Schweitzer et al. (1984) 
have shown by a recursive Green function method that in such a situation the density of 
states of the lowest magnetic subband agrees well with the exact result Eq. ( |3.7|) for the 
continuum model. It should be noted that this correspondence between the tight binding 
model and free electrons holds only for the lowest magnetic subbands near the edges of the 
spectrum. 
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D. Semiclassical Approximation and Network Model 



For potentials smooth on the scale of the magnetic length the localization in the integer 
QHE can be discussed in terms of semiclassical quantization and percolation.^ These ideas 
have been developed by Tsukada (1976), lordansky (1982), Kazarinov and Luryi (1982), 
Ono (1982a), Prange and Joynt (1982), Trugman (1983), Shapiro (1986), Wilkinson (1987), 
Chalker and Coddington (1988), Fertig (1988), Mil'nikov and Sokolov (1988), Jaeger (1991), 
and Lee, Wang, and Kivelson (1993). In the semiclassical limit, the motion of the electron 
separates into two components with vastly different time and length scales. On the one hand 
there is the slow motion of the guiding centers of the cyclotron orbits along equipotential 
lines of the disorder potential that is governed by classical drift equations. On the other 
hand there is the rapid cyclotron motion of frequency around the guiding center within 
the classical cyclotron radius = ^/2n + lie from the guiding center. The wavefunctions 
are non-zero on strips of width Rc around the path of the guiding center. They are thus 
located at the edges of the regions where V^(r) < E (see Fig. H). Near minima of the potential 
the wavefunctions encircle these minima, while near maxima of the potential they encircle 
the maxima. In both cases the orbits are closed and the states localized. Since V^(r) is a 
random potential the properties of the regions V{r) < E are described by the continuum 
model of percolation (Trugman, 1983) with the regions bounded by the equipotential lines 
corresponding to the clusters of the percolation model. At one energy between the minima 
and maxima the largest percolation cluster extends throughout the whole system and at this 
energy the wave functions are extended. For a disorder potential with a distribution that is 
symmetric this energy corresponds to the band center. 

The rapid cyclotron motion around the guiding centers and hence the extent of the 
wavefunctions perpendicular to the equipotential line become important whenever two orbits 
approach each other on a distance less than the cyclotron radius. This happens near saddle 



^Note that an expansion in is also an expansion in powers of h. 
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points of the potential. There tunnehng between the two orbits becomes important and 
effects of quantum interference can become observable. These effects are most important 
near the percolation threshold because there the percolating equipotential line is split up by 
saddle points into localized, closed contours. Thus while tunneling and quantum interference 
have a negligible effect in the tails of the density of states, they can in fact change the critical 
behavior of the system. 

Chalker and Coddington (1988) created a model that captures the physics of percolation, 
tunneling and quantum interference near the percolation threshold by introducing a network 
of saddle points connected by equipotential lines. Due to its quantum nature this network 
model shows critical behavior that is different from the classical percolation problem but 
coincides with that for short-ranged potentials (see section [VI I H ) . 

In the following we will first introduce the semiclassical picture of the QHE and then 
motivate the network model as a simplification of this picture. 

The semiclassical approximation is most conveniently derived by replacing the coordi- 
nates {x,y) by guiding center coordinates {X,Y) and relative coordinates {^,ri), given by 



x = X + (3.35a) 
y = Y + ri, (3.35b) 



and 



^ = --^{py-eAy), (3.36a) 
rj = -^{p^- eA^). (3.36b) 



From the equations of motion for the pure system {V{r) = 0), we get 



i = ^[Ho,^]=cu,rj, (3.37a) 
r] = UHo,v] = -^c^- (3.37b) 
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We see that rf) indeed rotate with angular frequency Uc around the guiding center. Due 
to the commutation relations of p and r both the guiding center coordinates and the relative 
coordinates obey canonical commutation relations, 

[i,v]=ill (3.38a) 
[X,Y] = -ill. (3.38b) 



With Eq. ( 3.36 ), the Hamiltonian Hq can be written in terms of ^ and rj, 



Ho = ^ {e + V') , (3.39) 

which expresses the degeneracy of the Landau levels as Hq does not depend on X and Y. 

In the presence of a disorder potential ^(r) the degeneracy is lifted. The equations of 
motion for the center coordinates are 

?■ P dV 

If the potential V"(r) is smooth on the scale Ic, we can replace V{x, y) by V{X, Y) and obtain 
a drift of the guiding center along equipotentials, 

f dV 

X = (3,41a) 
-If- 

In this limit, the eigenenergies are 

E = {n + l/2)huo^ + V{X,Y). (3.42) 

Due to the breaking of the time reversal symmetry by the magnetic field the guiding centers 
can move only in one direction determined by the magnetic field and the slope of the potential 
(Fig. I). 

The electron eigenf unctions can be approximated by 
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i,{u,v) = C{u)xn{v)e'^^''^''\ (3.43) 

where the variables u and v parametrize the distance along and perpendicular to the equipo- 
tential line, respectively, 

C\u) oc ^ , (3.44) 

and (f){u,v) is a gauge-dependent phase (Trugman, 1983). The semiclassical quantization 
condition is that the phase (j){u, v) must change by a multiple of 2tt upon going around a 
closed contour. This condition determines the allowed energies. From Eq. ( |3.43| ) it follows 
that the overlap between different states is exponentially small as long as the separation 
between the corresponding contours is large compared to the magnetic length. 

The problem of finding the wavefunction of largest extent is equivalent to the problem of 
finding the largest connected cluster in a continuum percolation problem (Stauffer, 1979). 
The "diameter" C,p{E) of this cluster grows on approaching the percolation threshold E,. as 

ep(^)cx|E-E,|-^^ (3.45) 

where the percolation critical exponent up is 4/3 (den Nijs, 1979; Black and Emery, 1981). 

In addition to the semiclassical drift of the guiding center coordinates tunneling has to 
be considered near saddle points of the potential, where two orbits approach each other. 
In order to study the influence of tunneling on percolation Chalker and Coddington(1988) 
introduced a network model. To motivate this model, let us look at an electron moving 
along an equipotential line far away, compared to Ic, from other trajectories. Then its 
wavefunction is given by Eq. ( p.43|) . This wavefunction also contains information about 
the degrees of freedom perpendicular to the trajectory that we are not interested in. In its 
place, we consider a simpler quantity, the complex function Z{u) of the coordinate u along 
the equipotential line, defined by 

arg(Z(M)) = arg(^(M, v = 0)) (3.46) 

and 
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\Z{u) 



dvij*{u,v)juip{u,v), 



(3.47) 



where ju is the M-component of the current density operator and u increases in the direction 
of net current flow so that the integral is positive. Since no current can escape perpendicular 
to the equipotential line Z{u) can only change its phase along the trajectory. This phase 
change depends on the gauge and arc-length of the trajectory in units of Ic- 

Let us now consider the region near the saddle points (Fig. There two degenerate 
wavefunctions meet and, due to tunneling, the description in terms of the Z{u) breaks down. 
Let us denote by Zi, Z2, Z^, and Z4 the values of Z{u) on the four arms some distance large 
compared to Ic away from the saddle point, where the Z{u) are still fair descriptions of the 
exact eigenfunctions. Due to the directedness of the trajectories two of these describe ingoing 
currents (say, Zi and Z2) and two describe outgoing (Z3 and Z4). Solving the Schrodinger 
equation near the saddle point gives a relation between the ingoing and outgoing currents, 

'zA 



'Z{^ 



M\ 



Current conservation constrains the 2x2 matrix M to unitarity \Zi 
This implies 



? + \Z2? 



(3.48) 



J = mUm, 



(3.49) 



where 



J 



(3.50) 



which has the general solution 

/ 

M = 



e*'^2 



^ ^ cosh 6 sinh 6^ ^ 



\ 



J 



e*<^4 



(3.51) 



sinh 6' cosh^^ 

with real 9, (p. The phases (pi can be set to zero by a suitable choice of gauge so that the 
saddle point is completely characterized by a real parameter 9. For 9 the saddle point 
behaves classically and the electron almost always goes from 1 to 4: \Zi\'^ \Z^\'^. For ^ 1 
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it also behaves classically, but now 1 goes to 3: \Z\\^ ^ l-^sP- In between the saddle point 
behaves quantum mechanically. For a special value 6^ the saddle point will be symmetric. 
At this value the probability to scatter from 1 to 3 and to 4 are the same. 6c is given by 

sinh^c = l, (3.52) 

which has the solution 9c, = ln(l + v^) ^ 0.8814. Fertig (1988) showed how 9 depends on 
the energy, 

sinh(^) = exp(-7r7/2), (3.53) 

where 

Vi is the saddle point potential, and Va is the second derivative of the potential in suitable 
coordinates. 

The network model consists of a square lattice of saddle points, each described by a ^ 
parameter (see Fig. |TOD. Two neighboring saddle points are connected by directed links 
representing the equipotentials such that at every saddle point there are two incoming and 
to outgoing links. Each link is characterized by the phase change of Z{u) along the link. 
Since the lengths of the links and hence the phase change depends on the random potential 
these phases are taken to be uniformly distributed random variables. In the original paper 
(Chalker and Coddington, 1988) the 9 parameter was fixed for all nodes in the network. 
Generalizations of the model where the 6''s were random variables were studied by Chalker 
and Eastmond (1993) and Lee et al. (1993). 

In contrast to the other models of disorder the network model does not deal with wave- 
functions but with transmission probability amplitudes. Thus it lends itself naturely to the 
study of properties of the transmission matrix through the system but it is not obvious how 
to get information about the multifractal properties of the system discussed in section [IX C 



Mil'nikov and Sokolov (1988) employed the semiclassical approximation to calculate the 
inverse localization length ^{E)~^ in an intermediate energy range, ro-(/c/o")^ <^ E T„. In 
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this energy range the equipotential hnes differ significantly from the percolating equipoten- 
tial line only near the saddle points, where the percolating equipotential line is split up. Due 
to the first strong inequality the parameter 7 of Eq. ( |3.54| ) is large at these saddle points 
and with Eqs. ( |3.48| ) and ( |3.53| ) the transmission probability through the saddle point is 
exponentially small. Mil'nikov and Sokolov assume that in this limit a single path, which 
consists of saddle points connected by parts of equipotential lines, is important for the trans- 
mission probability between points far apart in space. Then they treat this problem using 
the semiclassical theory of tunneling through a one-dimensional barrier. The solution to the 
eigenvalue equation (|3.42|) defines a dispersion relation Y{E,X) similar to the dispersion 
relation of a one-dimensional particle in a potential k{E,x). The solutions Y{E,X) are 
real on the equipotential lines but are complex in the classically forbidden regions near the 
saddle points. With Eq. ( p. 421) the inverse localization length ( |6.2| ) can be written as 



^E)-' = min|Imr„(E)|//2, (3.55) 



where Yn{E) are the roots of Eq. ( p. 421 ). The average distance between saddle points is given 



by the diameter C,p{E) ( |3.45| ). Thus we can rewrite Eq. ( p. 55] ), 



i{E)-'=ip{E)'^j dX\\mY{E,X)\/ll (3.56) 

where the integral runs over the classically forbidden range of X near a saddle point. Since 
both the width of this region and the maximum of |Imy(i?,X)| are proportional to E^/"^^ 
the integral in Eq. ( |3.56|) is proportional to E and with Eq. ( |3.45| ) we find 

i{E)-^ oc E'^^+i = E''/\ (3.57) 



While the result = 7/3 agrees with the numerical value ( |8.4| ) it is not clear whether it 
describes the same physical situation. The derivation of the result is based on the assumption 
that only a single trajectory between two points is important. While this assumption seems 
justified ii E ^ Tfj{lc/(yY and each saddle point transmits predominantly into a single link, 
closer to the percolation threshold the probabilities for transmission into both possible links 
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become comparable. In this energy range the two-dimensional character of the percolating 
equipotential line has to be taken into account. It is this limit that is described by the 
network model, which has (see Sec. |VIII A] ) the same scaling behavior as the random Landau 
matrix model. 



IV. FINITE-SIZE SCALING 

In this section we want to review the finite-size scaling theory as it was developed for 
thermodynamic phase transitions. Since the transitions between quantum Hall plateaus are 
believed to be zero-temperature quantum phase transitions we want to analyze the numerical 
data in section [VII1| in terms of this finite-size scaling theory. We base our discussion on the 



review of finite-size scaling by Barber (1983). 

For infinite systems, thermodynamic quantities can diverge at phase transitions. The 
"usual" scaling theory describes the behavior of the thermodynamic quantities near the 
transition in terms of scaling laws. In finite systems there are no phase transitions and no 
singularities in thermodynamic quantities. Nevertheless in the simultaneous limit that the 
temperature approaches the critical temperature and the system size goes to infinity the 
scaling laws of the infinite system are refiected in scaling laws for the finite system involving 
the system size. It is this finite-size scaling theory that we want to focus on in this section. 



A. Single-Parameter Scaling 

Let us consider a (i-dimensional thermodynamic system that is finite in at least one 
direction. Since we are interested in the QHE in two dimensions, let us assume d = 2. Then 
the system could have the geometry of a quasi-one-dimensional strip or cylinder, depending 
on the boundary conditions in the transverse direction, or it could be finite in both directions. 
Let us further assume that the infinite d-dimensional system has a phase transition but that 
the system with the finite length scale has no transition. An example of such a system would 
be the Ising model. 
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Let Poo{T) be a thermodynamic quantity of the infinite system that diverges at the 
critical point T^., 

P^{T)<xt-P, (4.1) 

as t = (T — Tc)/Tc 0. For the d = 2 Ising model this could be the susceptibility and the 
parameter T is the temperature. For a finite system this divergence will be rounded and the 
function Pl{T) of the finite system will have a maximum at a position Tm{L) that will in 
general be different from but will approach Tc as L —> oo. The approach is characterized 
by a shift exponent A defined by 



Tr 



oc L-^, L ^ oo, (4.2) 



where L is measured in multiples of some microscopic length scale of the system. 

The central idea of finite-size scaling is that the magnitude of the finite-size effects are 
determined by a single length scale i{T) that diverges at the critical point, 

e(T) cx r^ t 0. (4.3) 

This length scale is identified with the correlation length of order parameter fiuctuations. 
Thermodynamic quantities will depend on the ratio 

y = L/^{T). (4.4) 

For ?/ ^ 1, the system does not feel the finite size and intensive quantities will not depend 
on L. For ?/ ^ 1, the finite size rounds off the thermodynamic singularities associated with 
the phase transition. Thus y = 1 defines another characteristic temperature T*{L) of the 
system where finite-size effects become important. From Eqs. ( |4.3|j4.4| ) it follows that 



, L^oo. (4.5) 

If ^ (T) is the only relevant length scale the shift exponent A must also be given hj X = l/u 
(Barber, 1983). In the following we will not consider modifications to the finite-size scaling 
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analysis necessary to deal with the shift of the "critical" point in finite systems. The systems 
that we want to analyze exhibit a symmetry such that the maxima of scaling functions are 
always located at the critical point of the infinite system. 

The finite-size scaling hypothesis states that for finite L and T near Tc, 

Pl{T) = L^Qiy). (4.6) 

The exponent uj is determined by requiring that Eq. ( [4.6| ) reproduces Eq. ( [4.1|) in the limit 
L ^ oo. This requires 

Q{y) oc y-", y -> oo, (4.7) 

and 

UJ = p/u. (4.8) 
For y — i> at fixed L, Pl{T) has to go to a finite value and hence 

Q{y) const., y ^ 0. (4.9) 

For the particular case of the correlation length C,l of the finite system the finite-size 



scaling relation (O) becomes 

UiT) = LQ^{y), (4.10) 



since comparing Eqs. (|4.1| ) and (^^ ) yields p = u and hence u = 1. 

The phenomenological finite-size scaling theory for a particular thermodynamic quan- 
tity outlined above can be related to and derived from the scaling theory based on the 
renormalization group for infinite systems (Suzuki, 1977; Brezin, 1982; Barber, 1983). The 
renormalization group establishes a relation between the coupling constants K of a Hamil- 
tonian H and the coupling constants K' of a Hamiltonian H', in which all lengths have been 
rescaled by a factor b relative to H, 

K' = R(K), (4.11) 
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where R is an analytic function of K. At a phase transition there exists a non-trivial fixed 
point K* of the recursion Eq. ( |4.11| ), K* = R(K*). Near this fixed point the transformations 



( [4. lip can be linearized and the matrix dRa/dKp evaluated at K* has the eigenvalues 

K = by'^. (4.12) 

Using these eigenvalues Eq. (|4.11| ) can be rewritten as 

C = &^"a, (4.13) 

where C,a is a non-linear scaling field that is a regular function of the coupling constants K. 
In terms of the scaling fields the fixed point K' = K is given by = 0. 

The singular part of the free energy per degree of freedom Z^, responsible for the diver- 
gence of thermodynamic quantities at the phase transition, can then be expressed in terms 
of the scaling fields Ca instead of the coupling constants K. /^(Ci, • • • ; L) has the scaling 
form 

/,(Ci,...;L) = L-'^F(L^'^Ci,...), (4.14) 

where L is the length scale of the system that is reached after iterating the renormalization 
group equation (^4.11|) a number of times starting from the microscopic Hamiltonian, and 
F is an analytic function of its arguments. From the free energy (|4.14|) the scaling behav- 



ior of other thermodynamic functions can be obtained by differentiating the free energy 
appropriately, 

p,.(Ci,...) = L-g(L^^Ci,...). (4.15) 

This expression is equivalent, in the thermodynamic limit, to the single-parameter finite-size 
scaling hypothesis ([4.6|) if 



C, = Q^t + 0{f) (4.16) 

and 
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yi = i/^. (4.17) 

where we have used Eq. ( ^4.3| ). Furthermore, all other scaling fields Ci? > 2, have to be 
irrelevant, i.e., yi < 0, i > 2. We can interpret this statement in the following way. The 
Hamiltonian if of a system depends in general on a large number of coupling constants 
Ki. In the representation in terms of the scaling field (a it depends on a large number of 
scaling fields. However, the thermodynamic properties in the limit L — oo depend only on 
a small number of scaling fields (a that have positive scaling indices ya- In conventional 
phase transitions these are related to the temperature t ( [4.16|) and the symmetry breaking 
field h, 

ChOch + 0{ht'^). (4.18) 

If a coupling constant K shows single-parameter finite-size scaling it means that the 
coupling constant K{bL) of the rescaled system is a function of the coupling constant K{L) 
of the original system and the scale factor b only, 

K{hL) = f{h,K{L)). (4.19) 

In the continuum limit, 6 — >■ 1, this property allows one to define a (3 function that is a 
function of K only, 

mm^'-^. (4.20) 



The fixed points of the renormalization group are given by the zeros of the j3 function ( [4.201) . 
Near the fixed point the j3 function can be linearized, 

l3 = yKHK/K*), (4.21) 

where Eq. ( [4.12| ) was used. For > {yx < 0) the fixed point is repulsive (attractive). 
Comparison with Eq. ( [4.17| ) shows that the localization length exponent v is given by the 



inverse slope of the /? function at the critical point. 
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B. Corrections to Scaling 



The scaling fields (a with negative scaling indices Ha < are irrelevant in the thermody- 
namic limit L oo. However, for finite systems they lead to corrections to the asymptotic 
scaling laws and it can be necessary to take these into account when analyzing numerical or 
experimental data for finite systems. For two scaling fields, a relevant field Ci with yi > 
and an irrelevant field (2 with 1/2 < 0, the renormalization flow of the "coupling constants" 



Xi = L^^Ci and X2 = Ly^C,2 under increase of the length scale L is depicted in Fig. |TI|. The 
arrows indicate the direction of change of Xi when L is increased. The point xi = X2 = Q 
is the fixed point of the renormalization group. The flow corresponds to the down-hill flow 
near a saddle point in a potential. 

Near the fixed point a thermodynamic quantity Pl(Ci,C25 • • •) can be expanded in the 
arguments Xa- From Eq. ( |4.15|) we get 

Pi(Ci, C2, . . .) = L-{Qo + aiL^^Ci + a2L-\y^K2 + •.•)• (4-22) 

If the largest irrelevant scaling index y2 is sufficiently large compared to the next leading 
index it might be sufficient to take only these two indices into account and one arrives at 
a two-parameter scaling theory. If there are other scaling indices close to y2 it might be 
necessary to treat more scaling fields. 

For the saddle point fixed point shown in Fig. |TT] the flow along the line C2 = away 
from the fixed point is described by the single-parameter scaling relation 

Pi(Ci,0) = L-Q(L^^Ci,0), (4.23) 

where yi = 1/v. The flow along the line Ci = towards the fixed point is also described by 
a single-parameter scaling relation 

PL{QX2) = L-Q{Q,L'\y^\C2). (4.24) 



Eq. ( [4. 241 ) shows that the corrections to scaling due to the irrelevant scaling field C2 exhibit 



finite-size scaling at the fixed point Ci = 0. 
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C. Dynamic Scaling 



So far we have focused our discussion on the static properties of scahng quantities. We 



will now consider quantities varying in space and time. Analogous to Eq. ( |4.1| ), we expect 
in the infinite system the wave vector q and frequency uj dependence 

Poo(T,g,^) ocr^Q(ge(T),^ro(T)). (4.25) 

The relaxation time to(T) describes the rate at which the system relaxes towards its equi- 
librium state. It diverges at the transition and the dynamical exponent z is defined by 

ro(T) oc r oc t-''\ (4.26) 

In a finite system fluctuation on long length-scales are cut off by the system size L and from 
Eq. ( [4.25|) we get the dynamic finite-size scaling relation 



Pl{T, u) = L-Q {L/i{T),ujT,{T)) , (4.27a) 
= L^g (L^^^t, uj-^/'^h) . (4.27b) 



In the limit — we recover the static scaling relation Eq. 



V. SCALING BEHAVIOR AND FIELD THEORY 

In this section we discuss field-theoretical approaches to the localization problem and 
their relation to scaling theory. While field theory in the context of the integer QHE provides 
a framework for a scaling analysis, it cannot at present provide quantitative results suitable 
for a comparison with experiments. We will therefore restrict ourselves to an outline of the 
motivation and basic ideas. Thus this section lacks the technical details found in the other 
sections. 
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A. Anderson Localization 



We will start with a brief survey of the problem of Anderson localization in the absence 
of strong magnetic fields. For a more detailed review the reader is referred to Kramer and 
MacKinnon (1993). 

Using prior work of Thouless and co-workers (Edwards and Thouless, 1972; Thouless, 
1974; Licciardello and Thouless, 1975) and Wegner (1976), Abrahams et al. (1979) argued 
that the average dimensionless conductance g = Gh/e^ should obey a single parameter 
scaling relation. They calculated the (3 function of the conductance, 

in the limits g 1 and g ^ 1. In the localized regime, g <^ 1, the conductance decreases 
exponentially with system size, g = gQexp[—aL) and hence 

]imf3ig)<x\ng. (5.2) 

For large conductance, macroscopic transport theory expresses for a d-dimensional "cube" 
of side L the conductance G in terms of the conductivity a, 

G{L) = aL''-^ (5.3) 

so that 

\im^P{g) = d-2. (5.4) 

Assuming a continuous /3 function Abrahams et al. arrived at the flow diagram depicted 



in Fig. |T2[ The /? function has a zero only in d > 2 dimensions. They argued that the (3 
function has no zero in d = 2 dimensions since the first correction to /5 in l/^f is negative in 
perturbation theory (Langer and Neal, 1966; Abrahams et al, 1979). 

The same behavior is obtained when the disordered electron problem is mapped onto 
a non-linear a model (Wegner, 1979; Efetov et al, 1980). We will briefly sketch how this 
model is derived. We start by representing Green functions, like 
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G+{r,r',E) = {r\{z-H)-y), (5.5) 

z = E + ie, e > 0, as path integrals 

/ L)[$+]<i)+(r)$+(r') exp ( f £r^+{r){z - H)^+{r)) 

^' = 1^^) 

The sign of e ensures the convergence of the integral. To represent advanced as well as 
retarded Green functions we need to introduce two additional fields $~ and $^ with the 
opposite sign of e. Equation (|5.6|) can be interpreted as the statistical average (<l>+(r)$+(r')) 



with the partition function 

Z = y"D[$+]L)[$+]exp (^y"rfV$+(r)(2-/7)$+(r)^ . (5.7) 

Higher-order Green functions are given by higher-order correlation functions of the fields. 

In order to perform the average over disorder one needs to get rid of the determinant in 
Eq. ( |5.6D that depends on the disorder potential. One way of achieving this goal is the use 
of the replica trick (Wegner, 1979; Schafer and Wegner, 1980; Efetov et ai, 1980). Here the 
fields are replicated m times so that the denominator becomes [det(z — if)]™. In the formal 
limit m — s> 0, performed at the end of the calculation, the determinant becomes unity. 
Instead, the supersymmetry method (Efetov, 1983) can be employed where commuting 
as well as ant i- commuting (Grassmann) fields are used. Then the determinants for the 
commuting and anti-commuting fields cancel exactly. 

Using one of these schemes, the average over disorder in Eq. ( |5.7|) can readily be per- 
formed for the Gaussian white-noise potential (p.5|). Since the disorder potential couples 
two fields this leads to a term containing a product of four fields in the exponential. These 
quartic terms can be decoupled by a Hubbard-Stratonovic transformation to matrix fields 
Q{r). The effective Lagrangian for this matrix is (Pruisken and Schafer, 1981) 

1 



L[Q] = J d\ 



jTTQ'{r) + TT\n{E-Ho-Q) 



(5.8) 



where we dropped convergence-ensuring terms. This equation is the starting point of the 



derivation of the non-linear a model. It proceeds by evaluating the integral in Eq. ( |5.8|) in 
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saddle-point approximation, suitably decomposing the Q-matrix, and expanding the trace 
of the logarithm. The result is an effective Lagrangian 

1. 



L{Q) = J (fr 



--Tr(9g) 



(5.9) 



where the coupling constant 1/t is proportional to the conductance. 

The symmetry of the matrices Q depends on the symmetry of the original Hamiltonian 
and the choice of commuting or anti-commuting variables. 

For the non-linear a model ( |5.9| ) the (3 function for the coupling constant and hence for 
the conductivity can be calculated (Brezin et ai, 1980). The results depend on the symmetry 
class of the Q-fields. One distinguishes three universality classes: The orthogonal class for 
random potential scattering, the symplectic class in the presence of spin-orbit coupling and 
the unitary class in the presence of a magnetic field. The results of Wegner (1989) may be 
written in terms of the dimensionless conductance g as follows: 

(d - 2) - g-^ - |C(3)^~^ + 0{g-^), orthogonal case 
(d - 2) - 2g-^ - 6g~^ + 0{g-^), unitary case (5.10) 
(d - 2) + g"^ - f C(3)^~^ + 0{g-^), symplectic case 
Only in the symplectic case does the (3 function have a zero in d = 2 dimensions while in 
the other cases all states are localized. In particular, the presence of a magnetic field in the 
unitary case does not lead to extended states necessary for explaining the QHE. We will 
come back to this problem in the next section. 

Further support for the scaling picture of Abrahams et al. (1979) comes from numerical 
calculations (MacKinnon and Kramer, 1981; Pichard and Sarma, 1981; MacKinnon and 
Kramer, 1983). For a recent review see Kramer and MacKinnon (1993). Numerically it 
is easier to calculate the localization length instead of the conductance. This can be done 
recursively for strips or bars using either the Green-function method of section |V1 A| or 
the transfer-matrix method of section |V1B|. The numerical results are compatible with 



the finite-size scaling hypothesis ( [4.101 ). In one- dimensional systems no extended states 



are observed. In two-dimensional systems a transition is only observed in the presence of 
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spin-orbit scattering. In three-dimensional system a localization-delocalization transition is 
observed at a critical disorder. 

After the discovery of the absence of self-averaging in mesoscopic systems (Anderson 
et al, 1980; Altshuler, 1985; Lee and Stone, 1985) it became apparent that the description of 
the metal-insulator transition in terms of the average conductance is incomplete. Altshuler, 
Kravtsov, and Lerner (1986; 1989) argued that the whole distribution of the conductance 
needs to be considered. Single-parameter scaling in this context was discussed by Shapiro 
(1987). His scaling parameter need not necessarily be the average or typical conductance 
but any parameter that uniquely determines the distribution function. 



B. Nonlinear Sigma Model and Topological Term 

In the previous section it was pointed out that the absence of extended states in the non- 
linear a model for the unitary two-dimensional case apparently contradicts the explanations 
for the QHE (Aoki and Ando, 1981; Laughlin, 1981; Halperin, 1982) which require the 
existence of extended states for sufficiently small disorder. Pruisken (1984) and collaborators 
(Levine et al., 1983; Levine et al, 1984a; Levine et al, 1984b; Levine et ai, 1984c) argued 
that the failure of the non-linear a model is due to the emergence of a topological term in 
the effective Lagrangian that is not accessible to the perturbative treatment that leads to 



the Lagrangian (|5.9| ). A review of the discussion can be found in Pruisken's chapter in the 
book by Prange and Girvin (1987). Here we want to sketch the main features of this field 
theory. 

The Lagrangian (|5.9|) has to be augmented by a term proportional to the Hall conduc- 
tivity, 

C = Jd\ (-^o-LTr(9g)2 + ^a°,Trg [d.Q, dyQ]^ . (5.11) 



The values of the coupling constants cr°^ and a^y in ( |5.11| ) are the mean-field values corre- 



sponding to the conductivity tensor on short length scales. They are given by the SCBA 
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(Ando et ai, 1975). Upon increase of the length scale the conductivities get renormalized 
to their physical values on long length scales. We will use the correspondence between the 
scaling on short length scales and the SCBA in our discussion of the corrections to scaling 
in Sec. ( [VITrBD . 



The occurrence of two coupling constants 0"°^ and 0"°^ in the Lagrangian leads to a 
coupling of the renormalizations of axx and axy and changes the critical behavior compared 
to the non-linear a model for Anderson localization. In a dilute-instanton-gas approximation 
for the physically relevant field configurations in the weak-coupling limit axx ^ 1 Pruisken 
derived (3 functions for the renormalization with system size L of both cxxx and axy, 

d 1 

(3xx = = (T,.De-2— cos(27ra,,), (5.12a) 

Pxy = = -axxDe-^^^^=^ sm{27raxy). (5.12b) 

This leads to the renormalization-group flow diagram proposed by Khmel'nitskii (1983) 



shown in Fig. |T3|. It has two kinds of fixed points: stable fixed points at axx = 0, axy = n 
and fixed points at half-integer values of axy with a finite conductivity a^^. The former 
fixed points correspond to the localized wave functions of the model in the absence of the 
topological term while the latter correspond to the extended wave functions that carry the 
current in the QHE. Furthermore the flow diagram is periodic in axy so that the critical 
behavior is independent of Landau level index, since (Levine et ai, 1983) 

^ J d'rTiQ [dxQ, dyQ] = 2mqaly. (5.13) 

The delocalization fixed points at half-integer axy are partly attractive and partly repul- 
sive like the fixed point discussed in connection with Fig. [Tl|. They are characterized by a 
relevant scaling index yi = 1/v > Q associated with the deviation from 0"*^^ = {2n + l)/2 
and an irrelevant scaling index ?/2 < associated with the deviation from a*^. The starting 



point for the flow in Fig. 13 at microscopic length scales is given by the expressions for the 



conductivity tensor in the self-consistent Born approximation (SCBA) (Ando et ai, 1975). 
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While the field theory with the Lagrangian ( |5.11| ) provides a very appealing framework 
for the quantum Hall transition it has major shortcomings. Up to date it has not led 
to quantitative predictions, besides the existence of isolated critical fixed points. It has 
not been possible to calculate the longitudinal conductivity a*^ at the critical point, the 
localization length exponent z/, or the irrelevant exponent ?/2- The validity of the dilute 
instanton gas approximation that forms the basis of the approximate fiow equations (|5.12|) 
has been questioned by Weidenmiiller and Zirnbauer (1988). Recently, Zirnbauer (1994) has 



argued that the fiow-diagram Fig. |1^ does not describe the fiow of the coupling constants in 
the nonlinear a model ( p.ll| ) and that a different theory is called for. We conclude that at 
present the most promising tool to obtain quantitative results about the critical properties 
of the integer quantum Hall system seems to be numerical simulation, on which we will now 
focus our attention. 



VI. LOCALIZATION LENGTH 

The localization length ^{E) is one of the simplest quantities to calculate numerically that 
shows the occurrence of a metal-insulator transition. It can be defined by the asymptotic 
behavior of the modulus of the single-particle Green function 

G{r,r'-E) = {r\^^\r'). (6.1) 

The disorder average of G{r, r'; E) is a short-ranged function due to the average over the 
phase fiuctuations. It decays on the scale of the elastic mean-free path and does not 
show critical behavior. By averaging the modulus these phase cancellations are avoided 
and \G{r,r'] E)\ becomes long-ranged near the phase transition. |G(r, r'; ii^)| is not a self- 
averaging quantity and has a broad distribution. The inverse localization length 

^(E)-i = - lim , ^ ,A n\G{r,r';E)\ (6.2) 

r-r'|-»cxD |r — r I 

on the other hand is a self-averaging quantity so that it takes on the same value for almost 
all members of the disorder ensemble. The self-averaging property makes it particularly 
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useful for numerical calculations. In Eq. (|6.2| ) the energy E has to have an imaginary part 



since after averaging \G{r,r'; E)\ is non-analytic on the real axis. Using the self-averaging 
property and calculating the localization length for a single disorder realization in a finite 
system this imaginary part can be set to zero as long as E is not an eigenenergy of the 
system. 

In order to implement a finite-size scaling analysis of the localization length one studies 
long strips (in = 2 dimensions) or bars (in (i = 3 dimensions) and calculates the localization 
length along the infinite direction of the system. This localization length then depends on 
the perpendicular dimensions of the system (MacKinnon and Kramer, 1981; MacKinnon and 



Kramer, 1983). In section [VI A] we review a method to calculate recursively the localization 
length for systems of this geometry and extend it to the case of the random Landau matrix. 

Another way to define the localization length is to identify it with the inverse of the 
smallest Lyapunov exponent of the transfer matrix of the system (MacKinnon and Kramer, 
1981; MacKinnon and Kramer, 1983; Pichard and Sarma, 1981). We will return to this in 
section IVIB. 



A. Recursive Green Function Method 



The recursive Green function method was developed by MacKinnon and Kramer (1981; 
1983) for the Anderson model (p.31|) . In one dimension the corresponding Schrodinger 



equation (|3.32| ) reads 

Oj+i = {E - ei)ai - ai_i, (6.3) 

where the hopping matrix element was set to unity. Eq. ( |6.3D allows one to calculate re- 
cursively the wave function amplitudes aj. In higher dimensions the corresponding equation 
is 

A,+i = iE-H,)A-A_,. (6.4) 
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Hi is the Hamiltonian of the i^^ {d — 1) -dimensional shce of the d-dimensional "bar" and all 
quantities are {d — l)-dimensional matrices. Eq. ( |6.4| ) can be rewritten in terms of Green 
functions Gi^n coupling sites in slice 1 to sites in slice i, 

Gi,i+i = Ci.iCj+i.i+i, (6.5a) 
Gi+i.^+i = {E- - G,,,)-\ (6.5b) 



Iterating Eq. (6^) the matrix elements of the Green function connecting sites at both ends 
of long "bars" can be calculated without need for diagonalization of the Hamiltonian of the 
whole system. 



Eq. (6^) is tailored to the Anderson model with constant nearest-neighbor couplings. 



To apply this method to the Hamiltonian ( p.l6 ) it has to be modified. The random Landau 



matrix {nki\V\nk2) of Eq. (|A1|) describes a random one-dimensional tight-binding Hamilto- 
nian similar to the Anderson Hamiltonian but with long-ranged hopping matrix elements. 
To be able to apply recursion relations similar to Eqs. (|6.5| ) Aoki and Ando (1985; 1985) 
reinterpreted the one-dimensional Hamiltonian {nki\V\nk2) as the Hamiltonian of a two- 
dimensional system with couplings only between sites in neighboring slices. This is possi- 
ble since | (nA;i|\^|nfc2) P is essentially zero for \ki — k2\ > K and can be neglected, where 
is a constant of order {f3lc)~^ such that A''^^ = KLy/2TT is an integer. States with 
k E [{i — 1)K + l,iK[ belong to slice i and are coupled only to states in slices i — 1, i, 
and i + 1. {nki\V\nk2) can then be rewritten as a tri-diagonal matrix with diagonal ele- 
ments Hii being the Nk x Nk submatrices of {nki\V\nk2) connecting states within slice i 
and off-diagonal elements -ffj.i+i being the Nk x Nk submatrices connecting states within 
slice i with those in slice i + 1. The matrix elements in the off-diagonal matrices Hi i^i are 
now random numbers in contrast to the case of the two-dimensional Anderson model so that 



the recursion relation ( |6.5|) has to be modified to 



^1,1+1 — Gi^ji/ji+iGj+i^j+i, (6.6a) 
= (E — Hi^i — Hi^iiGi^iHi^i^i) (6.6b) 
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A different approach was introduced by Huclcestein (1990) who generahzed the one- 
dimensional Green function method to the case of long-ranged hopping matrix elements. 
Consider the random Landau matrix H^^^ of dimension K describing a system containing 
states k = 2ni/Ly, i = 1, . . . , K. In the following, we will label the states by integers 
k instead of quasi-momenta 2TTk/Ly. Adding the state K + 1 to the system leads to a 
Hamiltonian if (^+^) that can be decomposed into a block- diagonal part Hq and the couplings 
H' of state + 1 to the rest of the system, 

K+l 

ff{K+i) ^ J2 \k){k\V\k'){k'\, (6.7a) 

k,k'=l 

= Ho + H', (6.7b) 
Ho = H^''^ + \K){K\V\K){K\, (6.7c) 

K 

H' = J2i\k){k\V\K + 1){K + l\ + h.c.), (6.7d) 

k=l 

where we have suppressed the Landau level index. The related Green functions are 
{E - H^^^)-^ and Go = {E - Hq)-^ for which we obtain the recursion relation 



E-{K + i\v\K + iy (^-^"^ 

G^+i = Go + GoH'G^+\ (6.8b) 

The last relation can be used to calculate G^~^^ owing to the fact that the matrix H' only 
couples to site K+l. Using the property that the random Landau matrix is a banded matrix 
the dimensions of the matrices in Eqs. ( |6.8|) can be kept finite even in the limit K ^ (yo. 
Choosing matrices of dimension 2M/ /3, where M = Ly/ -\/27r/c, results in neglecting matrix 
elements of order exp(— 47r) compared to the diagonal matrix elements. Let be the 



dimension of these matrices, then the recursion relation (|6.8| ) for the matrix elements are 
(K + 1\G^^+^^\K +1) = 



K ^ -1 



E-{K+1\V\K+1)- Yl {K+l\V\i){t\G^'''>\j){j\V\K+l)] (6.9) 

i,j=K-N+2 

for the diagonal matrix element, 
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K 

(^IG^^+i^li^+l) = + {t\G^'^\3){3\V\K + l) (6.10) 

j=K-N+2 

for i < K, and 

= + E + + (6.11) 
for i < K and j < i^. 

The focalization length Aa/ of a system with width M = Ly/ \p2j^lc and length K is given 

by 

X-A^-.K) = --^ln|(l|GW(£;)|ir)|. (6.12) 

In calculating (l|GW|ir), we are faced with the difficulty that this is an exponen- 
tially decreasing quantity. To circumvent this problem we introduce variables Qk = 
{1\G^^^\K) / {1\G^^~^^\K -1) for the change in {1\G^^^\K) when the system size is increased 
by one state. In terms of the g^, Eq. (|6.12|) is given by 



c 

K 



(lIG^Iir) 


{1\G^^-^^\K- 


1) 


{1\G(^-^)\K- 


1) (1|G'(^-2)|K- 


2) 



■■■(1|G(^)|1) 



(6.13a) 



^ Eln|g,|, (6.13b) 



where we set qi = (1|G''^"'^^|1) = {E — (1|K|1))~"^. From Eq. ( |6.10D we get a recursion relation 
for qk+i, 

QK+i = (ir+l|G(^+^)|ir+l) E (l|^?W|j)(j|V|K + l), (6.14) 

j=K-N+2 

where (l|5'''^''|j) = (l|G'^^''|j)/(l|G'*''^^|-ft") is normalized to compensate for the exponential 
decay. With Eq. ( |6.11| ), we obtain the recursion relation for (l|5'^^'^^''|j), 

+ X: {l\g^''''\k){k\V\K+l){K + l\G^''+%)]. (6.15) 

k=K-N+2 / 
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B. Transfer-Matrix Method 



The transfer-matrix method deals with a physical situation where a long narrow disor- 
dered system is connected at both ends to semi-infinite ideal conductors. The transfer matrix 
T relates the amplitudes on the right-hand side of the system to those on the left-hand side. 
For the Anderson model (|6.4|) the transfer matrix Tj of slice i is given by 



E-H, -I 
I 



(6.16) 



where I is the unit matrix, since Eq. (|6.4| ) can be rewritten as 

A, \ 



T,. 



(6.17) 



For Chalker's and Coddington's network model of section [III D| the transfer matrix T consists 
of a product of four matrices, T = ABCD, where A and C are diagonal matrices with random 
phases, corresponding to the links connecting the nodes of the network, and B and D are 
block-diagonal, consisting of the 2x2 matrices M of Eq. ( |3.51| ). Note that the network of 



Fig. |T0| is periodic when using two columns of nodes as a unit. 

From the definition of the transfer matrix it follows immediately that the transfer matrix 
of a long system is just the product of the transfer matrices of the slices composing the 
system, 



1=1 



(6.18) 



where M„ is the transfer matrix of the system consisting of n slices and Tj is the transfer 
matrix of slice i. Oseledec's theorem (Oseledec, 1968) states that a limiting matrix F exists, 



hm (M;Mji/2n 



(6.19) 



with eigenvalues exp7i, . . . , exp72m, where 2m is the dimension of the transfer matrix. The 
7j are the characteristic Lyapunov exponents of M„. The inverse localization length is 
given by the Lyapunov exponent of smallest absolute magnitude, = min|7j|. 
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VII. OTHER SCALING QUANTITIES 



The localization length is by no means the only quantity that can be used in finite-size 
scahng studies. In fact, early numerical calculation on the QHE used the size-dependence of 
the Thouless number to obtain information about the critical behavior (Ando, 1983; Ando, 
1984a; Ando, 1984b). Information about the localization properties of eigenstates can also 
be obtained from their topological character. Huo and Bhatt (1992) used the density of 
states with non-zero Chern number as a finite-size scaling quantity. 

A. Thouless Number 

Edwards and Thouless (1972) have argued that the sensitivity of eigenenergies of finite 
systems to changes in the boundary conditions can be used to distinguish between extended 
and localized states. For locahzed states the eigenenergies should be insensitive against 
the change in boundary conditions if the system size is large compared to the localization 
length. Extended states should feel the change in the boundary conditions and the shift in 
energy should be of the order olh/r, where r is the time it takes the electron to diffuse to 
the boundary (Thouless, 1974). In second-order perturbation theory the change in energy 
5E when periodic boundary conditions are replaced by anti-periodic boundary conditions is 
approximately related to the conductivity a by (Licciardello and Thouless, 1975) 

where A£' is the mean level spacing at the Fermi energy and / is a constant of order unity 
depending on details of the microscopic model. g{L) is called the Thouless number. In 
the presence of a strong magnetic field Ando (1983) argued that the constant / is 7r/2. 
The Thouless number g{L) for a particular level is a strongly fiuctuating quantity as a 
function of level index. To get a less fiuctuating quantity typically the geometric mean 
5E = exp((ln5i5j)ave) is used, where the average is taken over a given energy interval. On 
the other hand, the very fact that g{L) is strongly fiuctuating was used by Fastenrath, 
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JanBen and Pook (1992) in a multifractal analysis of the integer quantum Hall transition 

(Sec. prui) . 



B. Chern Number 

While the Thouless number relates to the longitudinal conductivity a^x^ the ability of 
an eigenstate to contribute to the Hall conductivity Oxy serves as another means of distin- 
guishing extended from localized states. For a system with generalized periodic boundary 
conditions, 

t(L,)|V^)=exp(2^,)|V) (j = l,2), (7.2) 

where t(Lj) is the magnetic translation operator (Zak, 1964; Haldane and Rezayi, 1985; 
Arovas et al, 1988), the Hall conductivity may be written as a sum over all occupied states 
below the Fermi energy E (Thouless et al, 1982; Pook and Hajdu, 1987) 

,2 



X 



06. / \ dOo 



(7.3) 



= exp[-i{9ix/Lx + 92y/Ly)]\ij^), (7.4) 

where T* is the torus 6i G [0,27r[, z = 1,2. It has been shown by Thouless and co-workers 
(Thouless et al, 1982; Niu et al, 1985) and others (Avron et al, 1983; Kohmoto, 1985; 
Aoki and Ando, 1986; Pook and Hajdu, 1987) that the integral in Eq. ( |7.4| ) is 27r times 
an integer Ci(m), the first Chern index, which is a topological invariant characterizing the 
topological properties of the wave functions. Arovas et al. (1988) have shown that for 
states with finite Ci{m) the zeros of the wave function can be moved to any position in real 
space by a suitable choice of the boundary angles 61,62, while states with zero Ci{m) have 
zeros that are confined in space. This leads to a natural distinction between extended and 
localized wave function. Again, as in the discussion of the Thouless number, the sensitivity 
to boundary conditions is used to distinguish extended from localized states. Since the first 
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Chern index is an integer it allows a clear definition of extended versus localized states even 
for finite systems. 

VIII. NUMERICAL CALCULATIONS 

First numerical calculations of the critical properties for two-dimensional systems in a 
strong magnetic field were performed by Aoki (1977; 1978; 1985), Ando (1983; 1984a; 1984b), 
and Schweitzer, Kramer, and MacKinnon (1984). Mostly due to the insufficient computing 
power available at the time these studies were not able to obtain accurate information about 
the behavior of the localization length or the Thouless number. Ando (1983; 1984a; 1984b) 
interpreted his results obtained from the size dependence of the Thouless number as being 
compatible with an extended state at a single energy in each Landau level, a result previously 
obtained by Ono (1982b) by summing a certain class of diagrams analytically. However, the 
divergence of the localization length obtained numerically was much weaker than Ono's result 
^{E) oc exp(l/i?^), where the energy is measured relative to the critical energy. While Ando 
was not able to extract a result for the localization length exponent u from his data he found 
first signs of a single-parameter scaling law since his numerically constructed (3 function for 
the conductance depended only the conductance itself. Schweitzer, Kramer, and MacKinnon 
(1984) on the other hand interpreted their results for the localization length as indicating 
a band of extended states separated from localized states in the band tails by two mobility 
edges. 

In the following we will concentrate on more recent results that shed light on the critical 
behavior near the band centers of the Landau levels. In our discussion we have to keep 
in mind that numerical calculations always deal with finite systems while strictly speaking 
phase transitions exist only in infinite systems. Thus, we need a theoretical framework in 
which to analyze data for finite systems in order to extrapolate results to infinite system 
size. The finite-size scaling theory outlined in section |IV| provides such a framework. 

We will first review calculations for Hamiltonians projected onto a single Landau level. 
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In the strong magnetic field limit, fiujc S> F, this seems to be a reasonable approximation. 
The effects of Landau level mixing will be discussed in Sec. |V111D| while the influence of a 
periodic potential on a single Landau level is considered in Sec. |V111E| . 



A. Lowest Landau Level 

The first calculation that produced reasonable estimates of the critical exponents for 
the quantum Hall system was performed by Aoki and Ando (1985) (Ando and Aoki, 1985). 
They used a real-space model for the disorder with the random potential given by a sum over 
either short-range scatterers (|3.8|) or Gaussian scatterers ( |3.11|) . Using the recursive Green 
function method of section |VI A| they calculated the localization length for strips of width 



1.25 to 40.1 Ic- In the direction perpendicular to the strip periodic boundary conditions were 
employed to avoid the effects of edge states extended along the edges of the strip. Aoki 
and Ando tried to interpret their data in terms of Pruisken's two-parameter scaling theory 
(Sec. |V B| ) . However, in view of later calculations (Huckestein and Kramer, 1989; Huckestein 
and Kramer, 1990) it has to be realized that their systems were too small to be in the critical 
region. Nevertheless, they found that the localization length diverges at a single energy at 
the center of the Landau band and estimated the localization-length exponent to be z/ 2 
in the lowest Landau level. 

Huckestein and Kramer (1990) extended these calculations by introducing the random 
Landau- matrix model (Sec. |IIIB|) . Comparing their results for the localization length with 
the results of Aoki and Ando shows (Huckestein and Kramer, 1989) that the different random 
potentials give the same results for a width of 40 Ic, the largest calculated by Aoki and 
Ando, while for smaller widths differences are observed. Microscopic details of the disorder 
potential therefore seem to become irrelevant for the behavior of the localization length on 
this length scale for (5-correlated potentials. In Fig. ^ the normalized localization length 
\m{E)/M is plotted as a function of the system width M = Ly/ (27r)^/^Zc for energies between 
0.01 and 1.0 F relative to the band center. The length of the systems was 2.5 x 10^/c (2.5 x 10^/c 
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for the widest systems) and the width ranged up to 321 Ic. The data within the dotted region 
are compatible with a single-parameter scaling law 

Xm{E) = MA{M/^{E)), (8.1a) 
\E-E-'' 



aE) = eo 



(8.1b) 



At this point we want to elaborate on the significance of this statement. In Sec. [VI A| it 



was argued that the inverse localization length calculated from the Green function is self 
averaging for system length L —>■ oo. Numerically, it can only be calculated for finite L. 
The value obtained fiuctuates within the disorder ensemble. The statistical error can be 
estimated from sampling different realizations of the disorder. In addition, systematic errors 
of the order of C, {E) / L have to be expected. In the present calculation the statistical errors 
are typically less than 1%, with possible systematic errors of the same order. Using a 
statistical test (Huckestein, 1990), it was determined that the data in the dotted region are 
compatible with both scaling relations in Eq. ( p.l|) within the statistical errors. Corrections 
to Eq. (^ . 1|) larger than the statistical errors are ruled out by this procedure. The deviations 



from scaling behavior for smaller width in Fig. 14 are not only due to irrelevant scaling fields, 



as discussed in Sec. [IV B| , but are also associated with the particular one-dimensional limit 



of the system. As was pointed out by Ando and Aoki (1985) for M — the localization 
length can be calculated analytically. In this limit the distance AX = Akl"^ = 2Trl'^/Ly 
between the center coordinates of neighboring Landau states becomes large compared to Ic 
so that the matrix elements {nk\V\nk') fall off rapidly as a function of \k — k'\. The random 
Landau matrix becomes a tri-diagonal Hermitean matrix corresponding to a tight-binding 
model with nearest-neighbor hopping. Applying the exact result for the localization length 
of that model (Thouless, 1974) we get 

independent of energy. 

The critical region starts at system width M = 16, the same size at which the differences 
between the real-space model of Aoki and Ando and the random Landau-matrix model 
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vanish. For energies larger than O.SF no scahng is observed since the locahzation length 
does not decrease arbitrarily in the band tails but levels off at about one magnetic length. 



Fig. |T5| shows the data from the dotted region of Fig. |n| scaled according to Eq. (|8.1a| ) . Note 
that this figure contains data not only from Fig. |1^ but also from Figs. |17|, [1^ and from the 
network model (Sec. |V111 Ul ). The scaling relations Eq. ( |8.1| ) determine C,{E) only up to a 



constant prefactor. This prefactor can be fixed by observing that in the localized regime, 
M/^ ^ 1, Xm{E) converges to ^{E) for large M so that 

lim A(x) = x-^. (8.3) 

The divergence of the localization length in Fig. |1^ is given by the exponent 

1/ = 2.35 ±0.03. (8.4) 

This value for the exponent is in remarkable agreement with the experimental values 2. 3 ±0.1 
(Koch et ai, 1991b) and 2.4 ±0.1 (Wei et ai, 1994). By performing similar calculations for 
smaller systems of up to 150/c width Mieck (1990) obtained an exponent oi u = 2.3 ± 0.08. 

In the previous analysis the critical energy E^ was taken to be the band center so that the 
localization length diverges only at a single energy. The absence of extended states within 



a finite range of energies is consistent with the field-theoretic picture of Sec. [VB| and was 
proven for single Landau levels by Chalker (1987). 

The influence of a correlation length a = 1^ on the behavior of the localization length 
can be seen in Fig. |17| (Huckestein et ai, 1992). Single-parameter scaling behavior is again 
observed for the data in the dotted region. The necessary width increases by a factor of 
4 relative to the (5-correlated potential. One reason for this behavior is the increase in the 
effective length scale by a factor = (Z^ ±cr^)//^. Furthermore, due to the one-dimensional 
limit Eq. ( |8.2| ) the flow in Fig. |1^ starts at values of Xm/M smaller by a factor of 
compared to cr = 0. For large systems, M > 64, the scaling behavior is the same as for 
a = 0. Consequently, the data can be fitted to the same scaling curve in Fig. [T^ and the 



localization length exponent is again given by Eq. ( ^.4|) . The prefactor of the localization 
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length is roughly independent of the correlation length if the correlation-length dependence 
of the bandwidth = Tq/ (3 (Ando and Uemura, 1974) is taken into account. This result is 
in contrast to calculations by Ando (1989b) that showed a strong dependence of the prefactor 
on the correlation length. Presumably, the difference can be attributed to the limited system 
size M < 20 of that calculation. 

For significantly larger correlation length it seems difficult at present to reach the system 
width necessary to observe scaling behavior. For a = 4lc Mieck (1990) estimated a value for 
u of about 2.3. 

In another approach to the localization problem Huo and Bhatt (1992) calculated the 
Chern numbers of the eigenstates in finite systems using Eq. ( [7.4|) . By studying the size 
dependence of the disorder-averaged density of states with non-zero Chern numbers they 
found that a finite portion of the states in the band center is extended in the thermodynamic 
limit. Away from the band center they found that the width AE of the band of states with 
non-zero Chern number decreased with the area A of the system as 

AE oc A-^/^", (8.5) 

with u = 2.4 ± 0.1. Eq. ( |8.5| ) is obtained by assuming that the states are extended in an 
energy range where the localization length C,{E) oc \E — Ed^'^ exceeds the linear system size 
L, and A = L^. They obtain the same exponent from the scaling of average number of 
states with non-zero Chern number, oc A^, with y = 1 — l/2v. 

Ando (1992) studied the transmission through a disordered system with edges connected 
to two ideal leads. Using a lattice model for the system he calculated the transmission 
probability through the system as a function of the width of the system. He found that 
the energy range AE over which the transmission probability rises from zero to unity scales 
with the width M of the system as AE oc M^^^" with v = 2.2 ± 0.1. In the edge-state 
picture the finite conductivity near the center of the Landau levels is due to backscattering 
between edge states located at opposite sides of the sample. The backscattering is facilitated 
by extended bulk states and hence scales with the bulk critical exponent v. 
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Chalker and Daniell (1988) calculated numerically the two-particle spectral function 
S{q,uj] E) near the critical energy Ec- S{q,LJ] E^) is the Fourier transform of 



S{r, u- E,) = 6{E, - nuj/2 - E^)5{E, + nuj/2 - E^)7/>„(0)V'* (r)7/'^(r)^^(0). (8.6) 

Employing current conservation S{q,uj;E) can be expressed in terms of a wave-vector and 
frequency dependent diffusion coefficient D{q,u!), 

For a scale- invariant system, i.e., at the critical energy Ec-, the diffusion coefficient can 
only depend on the dimensionless quantity gL^, where = (hp{Ec)uj)~^^'^ is the size of a 
system with mean level spacing u. Comparing this result, which is supported by numerical 
calculations (Chalker and Daniell, 1988; Huckestein and Schweitzer, 1994), with the dynamic 
scaling relation ( |4.25| ) applied to the spectral function we can conclude that the dynamic 
critical exponent z = 2. 

In the limit q,uj the diffusion coefficient D{q,uj) becomes the diffusion constant Dq 
that is related to the conductivity by the Einstein relation 

= e'p{E,)Do. (8.8) 

Chalker (1987) has shown that due to a sum rule the diffusion coefficient D{qL^^) at large 
values of qL^ is reduced, 

D{qL^) cx (gL^)-^. (8.9) 

A non-zero value for rj is related to the multifractal properties of eigenfunctions at the 
mobility edge discussed in Sec. f^. Chalker and Daniell obtain values of rj = 0.38 ± 0.04 



and a^^ = 0.45e^//i. 

Huo, Hetzel and Bhatt (1993) improved the calculation of the conductivity at the criti- 
cal point using the same method as Chalker and Coddington. They found a universal value 
of a^^ = 0.5e^//i within uncertainties of about 10% for seven different kinds of disorder 
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potentials. Using the Kubo formula (Aoki and Ando, 1981) they calculated the Hall con- 
ductivity to obtain a*y = 0.5e^/h. While this result is trivial for random potentials, which 
are particle-hole symmetric on the average, they showed that it also holds for potentials 
with an asymmetric density of states, again with uncertainties of about 10%. 



B. Higher Landau Levels 

While the behavior of the localization length in the lowest Landau level seems to be 
well described by a single-parameter scaling relation the situation in higher Landau levels 
seems to be less clear. Published values for the localization-length exponent u in the second 
Landau level n = 1 range from 2.3 (Mieck, 1990) to 6.2 (Mieck, 1993). Indeed, it has been 
argued that the critical exponent u should be nonuniversal in higher Landau levels and 
should depend on details of the disorder potential, in particular the correlation length a of 
the disorder potential (Mieck, 1990; Mieck, 1993; Liu and Das Sarma, 1993; Liu and Das 
Sarma, 1994). Only recently have these discrepancies been reconciled with a universal value 
of u by considering corrections to single-parameter scaling behavior (Huckestein, 1994). 

The first estimate for the exponent u in the second Landau level was again given by 
Aoki and Ando (1985) who found u ^ 4. The disorder potential consisted of short-range 
scatterers in this case. The first systematic study of the effect of the correlation length 
of the disorder was performed by Mieck (1990; 1993). He noticed that the correlations of 
the matrix elements Eq. ( [AID in the random Landau model become independent of the 
Landau level index for large correlation length a (Huckestein, 1992). Mieck obtained values 
of 1/ ~ 6.2 for 0" = 0, z/ ~ 3.7 for cr = 0.5/c, and u ~ 2.3 for a = Ale- In performing a fit 



to the single-parameter scaling relation ( p.l|) for the short correlation lengths he used data 
from the energy interval between 0.5ro- and 0.8ra-. 

Huckestein (1992) showed that for a correlation length a = Ic the scaling behavior of 



the system is indeed the same as in the lowest Landau level. Fig. 18 shows the normalized 



localization length. The data from the dotted region can be fitted to the scaling curve 
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of Fig. O. Note that Fig. |T8| differs from Fig. O for the lowest Landau level only for 



small width. The behavior for short-ranged correlations with a = as shown in Fig. p 



is drastically different from both the behavior in the lowest Landau level (Fig. and the 



correlated potential (Fig. |18|). The important features of Fig. [T^ are the slow but noticeable 
decrease of \m/M as a function of M in the band center, the very weak energy dependence 
near the band center {E < O.SF), and the strong energy dependence towards the band tails 
{E > O.SF). The first observation shows that the system width is too small to reach the 
critical region of single-parameter scaling as \m{Ec) / M = Ac is constant at the critical 
point. This prohibits an analysis of the data in terms of a single-parameter scaling relation. 
If the strong energy dependence in the tails of the band is still analyzed using a single- 
parameter scaling relation large, nonuniversal values for u are obtained (Aoki and Ando, 
1985; Mieck, 1990; Huckestein, 1992; Mieck, 1993; Liu and Das Sarma, 1993; Liu and Das 
Sarma, 1994). In view of the absence of scale invariance in the band center and the fact 
that these energies are mostly outside the critical energy range in the lowest Landau level 
it seems highly doubtful that these exponents describe the critical behavior of the system. 
The weak system-size dependence in the band center is due to corrections to the single- 



parameter scaling laws as discussed in Sec. [IV B|. To analyze these corrections we have 



to connect the behavior in the band center for short correlation length to the universal 
single-parameter scaling at larger correlation length where the corrections are too small to 
be observed. From the dependence of Xm{Ec)/M on and M in Fig. ^ we see that the 



asymptotic single-parameter scaling regime is reached at about (3"^ = 1.8. The corrections 
for smaller values of show single-parameter scaling as a function of and M as expected 
from two-parameter finite-size scaling with one relevant and one irrelevant scaling index (see 
Sec. [IV Bp (Huckestein, 1994). Eq. ( [4.24[ ) for the localization length Xm reads 



Xm{Ec,0^) = MA (0,M-l^'-ICi„) , (8.10) 

where A{{M/^(E)Y^'^, 0) is the single-parameter scaling function of Eq. ( [8.1[ ). The physical 
nature of the irrelevant scaling field (i„ is unclear, except that it depends in some way on 
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the disorder potential, particularly on the correlation length parameter In the present 
context Eq. ( |8.1CI| ) is an ansatz that describes the numerical data well. 

Relating the unknown irrelevant scaling field Q„ to an irrelevant length scale ^irr oc 



This is a single-parameter scaling law similar to Eqs. ( ^.1|) but this time describing the 
fiow towards the critical point as indicated by the negative scaling index. By applying 



and the dependence of ^irr on the correlation length of the disorder potential. Fig. pl| 
shows \m/M as a function of M/^i„. The irrelevant length scale increases by four 
orders of magnitude when the correlation length of the disorder potential is reduced from 
0.8/c to (Fig. ^2]). The overall scale of is not given by the fitting procedure. It is 
approximately known by observing that is becomes of the order of the magnetic length 1^ 
as the correlation length approaches Ic- The asymptotic single-parameter scaling governed 
by the localization length ^ can only be observed for length scales larger than the irrelevant 
length scale ^irr. Systems numerically accessible at present are 3 to 4 orders of magnitude 
too narrow to reach the asymptotic single-parameter scaling regime. The effects seen in all 
previous numerical calculations for short-ranged potentials in the second Landau level do 
not reflect the asymptotic single-parameter scaling regime but are dominated by corrections 
due to irrelevant scaling flelds. The fact that these corrections can be analyzed using a two- 
parameter scaling relation shows that the system sizes accessible in numerical simulations 
are in fact sufficient to reach the critical regime. It is just the description in terms of 
single-parameter scaling laws that becomes valid only after the irrelevant length scale is 
exceeded. 

For small values of M~'^'"'(^ii.r we can expand the function A in Eq. ( p.lOp to get 



we can rewrite Eq. ( |8.1CI| ) as 





the same fitting procedure as in Sec. |VIII A| we can obtain the scaling function A(0, x) 



A = Ae + 6M-'j''"lCirr + ---. 
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From the width dependence of A — Ac in Fig. we get the irrelevant scahng index y-^^^ = 
—0.38 ± 0.04 (Huckestein, 1994). For M/^i„ ^0.1 this exponent describes the data well. 
This suggests that, at least in this regime, corrections due to further irrelevant scaling fields 
with scaling indices less than yi„ are negligible. 

The structure of the fixed point described by the two-parameter scaling relation 

Am(AE, Cirr) = MA (^M'/'AE, M-l^^"ICirr) (8.13) 

is identical to the one expected by Pruisken (1987) for the delocalization fixed points at 
half-integer axy (see Fig. |l^), 

AM(Aor,„ Ad,,) = MA (Mi/'^Aa,,, M-l^lAa,,) . (8.14) 

Here Aaij are the linear deviations from the fixed-point values of cxjj. Near the critical 
point Aaxy is proportional to AE so that the scaling with respect to the first variables in 
Eqs. ( |8.13| ) and ( ^.14| ) is the same. For the second variables, (i„ and Aa^x a direct relation 
is not obvious. Though both quantities depend on the correlation length of the potential it 
is not clear that Cirr is proportional to Aaxx- 

We can understand the dependences of the corrections to single-parameter scaling on the 
correlation length of the disorder potential and the Landau level index in the field-theoretical 
framework. As discussed in Sec. |VB| the mean-field approximation to the field theory is given 
by the non-critical SCBA (Pruisken, 1987). To observe the localization-delocalization fixed 
point brought about by fluctuations around the mean-fleld saddle point the system size has 
to exceed the localization length in SCBA. The scaling in this approximation is described 
by the corresponding unitary non-linear cr-model without topological term and is given by 
the (3 function in Eq. ( |5.10| ). The localization length is obtained by integrating the leading 
term in Eq. ( p.lO[ ) for the unitary case. In terms of the conductivity a^^ on small length 
scales it is given by (Wei et al, 1985) 

eo = /cexp(7rVL2). (8.15) 
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Using the result for the conductivity in SCBA for short-range scatterers (Ando and Uemura, 
1974) the dependence on Landau level index n is 

io{n) = /cexp [(2n + 1)^] . (8.16) 

While this length scale is of the order of the magnetic length in the lowest Landau level it is 
more than 3 orders of magnitude larger for n = 1. This explains why the corrections are not 
observed in the lowest Landau level even for 5-correlated potentials. The strong variation 
of the irrelevant length scale in Fig. can also be understood from the range dependence 



of the conductivity in SCBA and Eq. ( ^.15|) . For n = 1 the conductivity decreases by about 
a factor of 6 between a = and cr = l^. The corresponding length scale decreases from 
about 10^/c to about one magnetic length. In the lowest Landau level the conductivity 
changes by a factor of 2 and decreases only by about the same factor. 

Direct calculations of the critical value of the longitudinal conductivity axx are missing 
in higher Landau levels. Ando's results for the Thouless numbers for a lattice system with 
several Landau levels do not show a strong Landau level dependence as expected from SCBA 
(Ando, 1989a). It has been suggested that zero-temperature quantum phase transitions in 
two dimensions should have a universal value for the conductivity at the transition (Fisher 
et al, 1990; Fisher, 1990; Lee et ai, 1992). Lee, Wang, and Kivelson (1993) argued that the 
critical value a*^ is related to the fixed-point scaling amplitude Ac. They obtain a*^ = 1/2 
for Ae = 1/ ln(l + ^/2) which is very close to the best-fit result Ac = 1.14 ±0.02 (Huckestein, 



1994). Since the scaling curve Fig. |T^, and in particular Ac, is independent of Landau 
level index, this would support the notion of a universal value of a*^ = 1/2 in accordance 
with the numerical calculations in the lowest Landau level (Huo et al, 1993). Landau level 
independence of the critical conductivity is also a feature of Khmel'nitskii's flow diagram 
Fig. O. 
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C. Network model 



The first single-parameter scaling function for the QHE was obtained by Chalker and 
Coddington (1988). They studied the scaling behavior of the network model described in 
Sec. [1IID| using the transfer-matrix method of Sec. |V1 B| . Their scaling function is not 



reproduced here as it is identical to Fig. |15|. In fact, Fig. [T^ was obtained by jointly fitting 
data from the network model and the random Landau-matrix model for n = 0, a = 0, 
n = 1, a = Ic, and n = 1, a = Ic, to a single scaling function. The exponent u was estimated 
by Chalker and Coddington to be about 2.5 ± 0.5. The large uncertainties are due to the 
non- linear dependence of the parameter 6 on energy, Eq. (|3.53|) . From the fact that the 
scaling functions of the models are identical it follows that the exponents are identical, too. 

The network model has the remarkable feature that it shows very small finite-size cor- 
rections to single-parameter scaling. In contrast to the systems discussed above the data 
are close to the scaling function Fig. |15| already for small system width. This behavior can 
be understood from the different one- dimensional limits of the models. In this limit the 
normalized localization length of the network model is given by (Jaeger, 1991) 

^*^^^)-6n2 + lncoshf^)r\ (8.17) 



M V V 2 

where 7 is related to the energy E by Eq. (|3.54|) . In contrast to Eq. ( ^.2|) the network model 



retains an energy dependence in the one-dimensional limit and the value at the critical 
energy \m{Ec)/M = l/ln2 is only about 20% larger than the fixed point value Ac. 

A generalization of the network model in which the ^'s are taken to be random variables 
was discussed by Chalker and Eastmond (1993) and Lee, Wang, and Kivelson (1993). From 
finite-size scaling studies they find that this generalization is an irrelevant perturbation of 
the original model. If the distribution of the 6''s is narrow the scaling behavior is unchanged. 
Lee, Wang, and Kivelson obtain u = 2.43 ± 0.18 in this case. In the limit that the width 
of the distribution of 6''s goes to infinity the system behaves classically since each saddle 
point transmits only into one of the two possible links. Lee, Wang, and Kivelson find single- 
parameter scaling with u = 1.29 ± 0.05 compatible with the classical percolation exponent 
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4/3. The classical percolation fixed point is unstable against quantum effects and the system 
flows towards the quantum flxed point described hj u = 2.35±0.03. Close to this flxed point 
the corrections to single-parameter scaling due to the flnite width of the distribution of 6''s 
can be analyzed in the same way as discussed above. It was in this context that Chalker and 
Eastmond (1993) pioneered the analysis in terms of irrelevant scaling flelds in the context of 
the QHE. They obtained the irrelevant scaling index = —0.38 ± 0.02 in agreement with 
the results in the previous section. However, for the network model the corrections to the 



fixed-point value Ac are negative as opposed to Fig. 21 where they are positive. 



D. Effects of Landau Level Coupling and Spin-Orbit Interaction 

When the strength of the disorder is not negligible compared to the cyclotron energy huc 
one has to go beyond the single-Landau-level approximation and include the effect of Landau 
level coupling by the disorder potential. In this context two different questions are important. 
For weak Landau level coupling such that each Landau level still contains a critical energy, 
does the Landau level coupling change the critical behavior at these critical energies? On the 
other hand, for strong disorder the system is expected to become insulating eventually (cf. 
Sec. [V A|) . In this limit the question is: How do the divergences in the localization length 
disappear with increasing disorder? Numerically these question are hard to answer since 
the inclusion of several Landau levels necessitates much smaller systems so that finite-size 
effects become very prominent. It is also not obvious exactly how many Landau levels have 
to be taken into account to capture the essential physics. 

The first question has been approached by Liu and Das Sarma (1994). They considered 
the two lowest Landau levels and studied the scaling behavior of the localization length as 
a function of the strength of the Landau level coupling for short correlation length of the 
potential. They find that while in the lowest Landau level the critical exponent u does not 
change when the Landau level coupling is included and is close to the value of Eq. (p.4|), in the 
n = 1 Landau level there is a change in the exponent. While Liu and Das Sarma erroneously 
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concluded from these data that Landau level coupling is essential for universal scaling for 
short correlation length, we can understand this behavior by considering the influence of 
Landau level coupling on the magnitude of corrections to scaling. With increasing Landau 
level coupling the inter-Landau level scattering increases and the conductivity on short length 
scales decreases. From Eq. ( ^.151 ) we see that this drastically decreases the characteristic 
length scale for single-parameter scaling. The system thus moves closer to the quantum 
critical point and finite-size corrections become smaller. 

The dependence of the critical energies on Landau level mixing has been studied by Ando, 
both for continuous systems (1984b) and for lattice systems (1989a), using the Thouless 
number. For (5-function scatterers and a Hamiltonian projected onto the three lowest Landau 
levels he finds that the localization length diverges at critical energies in each Landau level 
and that these energies move upward in energy when the Landau level coupling is increased. 
It also has to be pointed out that when the level broadening becomes comparable to the 
Landau level separation the localization length becomes very large, even for energies between 
the original Landau levels. For lattice systems the critical energies also move upward and 
eventually vanish. The square-lattice system studied by Ando is symmetric with respect to 
the center of the tight-binding band and for every Landau level below this center there exists 
another Landau level above with the opposite Hall conductivity. When two Landau levels 
merge in the center of the tight-binding band their contribution to the Hall conductivity 
vanishes and so does the associated critical energy. Due to this symmetry for the lattice 
system it is not clear whether the upward motion of the critical energy with increasing 
Landau level coupling is a generic feature. 

The influence of spin degeneracy on the critical behavior has been studied by Lee and 
Chalker (1994) and Wang, Lee, and Wen (1994b) in an extension of the network model. 
To account for the spin degree of freedom two quantum-mechanical fluxes are carried by 
each link. The random phases on the links are then replaced by random U{2) matrices. 
The numerical results of both calculations support the existence of two separate plateau 
transitions even in the spin-degenerate case. The critical exponent of each transition is 
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in agreement with Eq. ( p.4|) for single Landau levels. If, on the other hand, the data are 
analyzed assuming a single critical energy over a limited range in energy a fit with an 
exponent v ^ 5.8 is possible. 



E. Presence of a Periodic Potential 



In this section we want to consider the infiuence of an additional periodic potential on 
the scaling behavior of the integer QHE. For definiteness, let us consider a periodic potential 
of the form 



In the absence of disorder each Landau level splits into p subbands if a = 27r/^/a^ = q/p (Az- 
bel', 1964; Hofstadter, 1976) and the Hall conductivity in the gaps between these subbands 
is quantized in integer multiples of e^/h (Thouless et ai, 1982). While the contribution of 
each subband to the Hall conductivity takes on large positive and negative values if p is 
large, the sum of these contributions over a whole Landau level equals unity. We see that 
at least in the absence of disorder the additional periodic potential is quite relevant in that 
it completely changes the phase diagram and produces phases with new quantized values of 
the Hall conductivity. Under the infiuence of disorder the subband gaps will eventually be 
filled in and the intricate structure of the phase diagram will evolve into the simple phase 
diagram of the integer QHE Fig. |l^. 

The phase diagram in Fig. ^ was obtained by performing finite-size scaling studies for a 
system with a = 3/5 as a function of the ratio T/Eq of the strength of disorder to periodic 
potential (Huckestein, 1993). The Hall conductivity was calculated from the Diophantine 
equation 



where a is the Hall conductivity in units of e^/h, t labels the subband gap, and s is an 
integer (Thouless et ai, 1982). In the presence of disorder the Hall conductivity can only 




(8.18) 



pa + qs = t, 



(8.19) 
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change when a critical energy moves through the Fermi energy. For strong disorder there 
is only one critical energy at the center of the Landau level where the Hall conductivity 
changes from zero to one, while for weak disorder all phases of the non-disordered system 
are recovered. 

The preceding discussion raises the question whether or not a periodic potential is a 
relevant perturbation of the fixed point at the center of the Landau level. For fixed system 
width the ratio Xm{Ec)/M increases with increasing strength of the periodic potential and 
decreases with increasing system width for fixed Eq (Huckestein, 1994). The analysis of 
this behavior shows that it is compatible with the scaling due to an irrelevant scaling field 



according to Eq. (|8.11|) . The scaling index obtained in this way is within error bars the 
same as in the second Landau level and in the network model, y^j-j. = —0.38 ± 0.04. Thus 
the periodic potential is irrelevant at the fixed point at the center of the Landau level. The 
universality of the irrelevant scaling index suggests that further scaling indices are sufficiently 
small compared to ?/irr so that they are not observed in the numerical simulations. 

The critical behavior at the other fixed points introduced by the periodic potential is 
still an open question. Numerically it is a hard problem since the finite-size corrections to 
single-parameter scaling become very strong and the position of the critical energy is not 
known a priori as it is for the critical state at the center of the Landau level. Since the new 
fixed points are connected to the band-center fixed point[| (see Fig. |2^) one could conjecture 
that the critical behavior is the same for all fixed points. 

IX. MULTIFRACTAL ANALYSIS 

So far our discussion focused on quantities, like the localization length, that are self- 
averaging in the thermodynamic limit, even at the localization-delocalization transition. 
These quantities are a rather special class. In general, phase transitions are accompanied by 



^At least for certain values of a. 
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large fluctuations in physical observables. At the transition where the natural length scale, 
the localization length, diverges these fluctuations appear on all length scales. In this section 
we want to analyze critical fluctuations of local observables in terms of multifractal measures. 
We will flrst outline the language of multifractal analysis, then apply this apparatus to 
computer simulations, and finally put this discussion into the context of the scaling theory. 



A. Generalized Dimensions and Singularity Strengths 

The multifractal analysis of measures was pioneered by Mandelbrot (1983) and further 
developed by, among others, Hentschel and Procaccia (1983), Halsey et al. (1986), and 
Chhabra and Jensen (1989). A recent review in the context of the localization-delocalization 
transition was written by JanBen (1994). 

Consider a normalized measure on a two-dimensional square of linear dimension L, 



Jd^ = l. (9.1) 

L2 



An example of such a measure is the local density, i.e., the squared modulus |'?A(r)p of an 
eigenfunction of the Hamiltonian. The large fluctuations of this quantity can be seen in 



Fig. ^ where the local density for a state at the center of the lowest Landau level is shown. 

With this measure we can deflne box-probabilities p{l) by integrating the measure over 
squares of linear dimension /, 

p{l) = J dfi= J rfV|^/'(r)|2. (9.2) 

n{i) Q{i) 

In our example p{l) is just the probability to flnd an electron in the box The dimension 
of the support of the wavefunction is ci = 2, since the wavefunction has only isolated zeros. 
By covering the system with squares we get A^(A) = A^'^ different box-probabilities 
Pi{l), where X = l/L. Due to the normalization Yl^Ji^ Pi{l) = 1 the average box-probability 
scales with the ratio A as 

(p(O)l oc A^ (9.3) 
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where the average is defined as 



= ^ g ■4.. (9.4) 



The scahng law (|9.3|) is not useful in distinguishing between localized, extended, and 
critical wavefunctions as it holds for all normalized wavef unctions. This situation changes if 
we look at the scaling of moments of the box-probabilities. In the absence of length scales 
these will also show power law scaling but with nontrivial exponents T(g), 

oc A"+^(^). (9.5) 

Here q is real, but not necessarily integer. For a homogeneous measure, Pi{l) = p{l) = 
1/N{X), one finds r(g) = {q — l)d. We therefore define generalized dimensions D{q), such 
that 

Tiq) = iq-l)Diq). (9.6) 



We need to clarify under which conditions we expect power law scaling as in Eq. ( p.5|) to 
hold. As already pointed out, the localization length ^ must not set a relevant length scale. 
For finite system size L this gives the condition L ^ ^. For short length scales microscopic 
lengths become important. In the quantum Hall system wavefunctions are smooth on scales 
of the magnetic length or the correlation length of the disorder. In order to see power 
law scaling we therefore need ^ where 1^ is a microscopic length. In summary, the 

condition for power law scaling is 

lm<^l < L (9.7) 

In the limit L oo this condition is fulfilled only at the critical energy where the localization 
length diverges. In this limit the generalized dimensions are given by 

{q-l)D{q) = liJ-^^^^4^-d. (9.8) 
A^o In A 

For finite system size the condition ( p.7| ) is met over a finite range of energies. In this case 
the generalized dimensions can be estimated from the slope of ln{p'^{X))L vs. In A over a 
range of parameters that satisfies Eq. (|9.7|) . 
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For g = we find from Eq. ( |9.5| ) tliat -D(O) = d is tlie dimension of the support of the 
measure. The exponent D{2) is related to the scaling exponent d* of the inverse participation 
ratio, defined by Wegner (1980) 

p(2) = [ rfV|^(r)|^ocL-'='*. (9.9) 



Formally, this equation is obtained from Eq. ( |9.5|) for g = 2 in the limit / — > 0. In this limit 
d* = D{2). In the limit g — 1, p^'^"^-* can be expanded in Eq. (|9.8|) and Di is given by 

A-o (p(A))LlnA 

D{1) is sometimes called the information dimension. 

The function D{q) has the following analytic properties that we will only sketch here. 
For a complete discussion see (Hentschel and Procaccia, 1983; Halsey et al, 1986). D{q) 
is a monotonously decreasing, positive function of q bounded from above and below by 
-D_oo = D{—oo) and D^o = -D(oo), respectively, the generalized dimensions that describe 
the scaling of the minima and maxima of the measure. 

Another way of expressing the multifractal properties of a measure is the spectrum of 
singularity strengths /(a). It is related to the function r(g) by a Legendre transform, i.e., 
we change from the variable q to the variable 

f{ot{q)) is then implicitly given by 

f{a{q)) = a{q)q-r{q). (9.12) 

From the analytical properties of the Legendre transform and the function r(g) one finds 
(Chhabra and Jensen, 1989) that f{a) is a positive function of negative curvature with a 
single maximum, a ranges from D^o to -D-oo- At these points /(«) vanishes with infinite 
slope. For g = 1, /(a) has slope 1 and /(a(l)) = a(l) = D{1). It takes on its maximum 
value /(ao) = for g = 0, where = a(0) > d. 
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Chhabra and Jensen (1989) introduced a method to calculate /(a) directly. It has the 
advantage to be numerically more stable than calculating r(g) and performing the Legendre 
transform ( |9.12| ). In addition to the box-probabilities we need to define a one-parameter 



family of normalized probabilities 

/i.(g,0 = [Pi(0]VE[Pi(Or, (9-13) 

i 

that act as g-microscopes on the original measure. For each value of q the corresponding 
values of a{q) and /(g) = f{(y{q)) can be calculated from 

a{q) = — (9.14) 

and 

Kq) = ^^i^. (9.15) 

Geometrically, the /(a) has a quite intuitive interpretation. If a certain box-probability 
Pi{l) scales with a singularity strength a, i.e., Pi{l) oc 1°", then /(a) is the fractal dimension 
of the set of boxes with singularity strength between a and a + da. 

Another aspect of the /(a) spectrum becomes clear when the average ( |9.4| ) is interpreted 
as an average over an ensemble instead of an average over the whole measure (Pook and 
JanBen, 1991; Fastenrath et ai, 1992; JanBen, 1994). In terms of the distribution function 
P{p, A) the average of a function A{p) can be written as 

{A{p,X))= f'dpP{p,X)A{p). (9.16) 
Jo 

Changing the integration variable to a = Inp/lnA and defining P{a,X)da = P{p,X)dp, 
Eq. ( |9.16| ) becomes 



{A{p,X))= daP{a,X)P{X''). (9.17) 
Jo 

From this we recover Eq. ( |9.5| ) for the scaling of the moments of p if (Pook and JanBen, 
1991) 
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P(a, A) oc A-^(")+^ (9.18) 

where /(a) is given by Eq. ( |9.12| ). We see that /(a) describes the scahng of the whole 
distribution of box-probabihties, while the r(g) describe the scaling of certain moments of 
that distribution. 

The connection between the shape of the /(a) curve and the distribution function P{a, A) 
becomes more transparent if we approximate the maximum of the f{a) spectrum by a 
parabola, 

Note that this parabolic approximation depends only on one parameter ao, besides the 
dimension d of the support of the measure, which follows from the property /(a(l)) = a(l). 
The corresponding probability distribution of the box-probabilities is log-normal. 

We see that the absence of length scales that is reflected in the /(a) spectrum corresponds 
to broad distributions of local quantities like the box-probabilities. The parameter that 
defines the parabolic approximation is the value of a at the maximum of /(a) and describes 
the scaling of the typical value ptyp(A) of the box- probability defined by 

Ptyp(A) = exp((lnp(A))) oc A"°. (9.21) 

The multifractal density fluctuations lead to anomalous diffusive behavior (Chalker and 
Daniell, 1988). This is reflected in power-law decay of density correlations and slow decay 
of temporal wave-packet autocorrelations (Huckestein and Schweitzer, 1994). The exponent 
7] of density correlations. 



|V'(r)^(r + R)|2 oc |R^^ (9.22) 
is related to the generalized dimension D{2). If we replace in Eq. ( |9.4| ) the average by an 



ensemble average and express p^(/) as an integral over (|9.22| ), we get 
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p^{l) (X (9.23) 

Comparing with Eqs. ( p75|) and (|9.6| ) yields the result 

ri = d-D{2). (9.24) 

At finite wave vector and frequency the exponent rj describes the decay of the wave-vector 
and frequency- dependent diffusion coefficient (Eq. (|8.9|) ). The value rj = 0.38±0.04 obtained 
by Chalker and Daniell (1988) from this quantity agrees well with Eq. ( p.24|) and D{2) = 
1.62 ± 0.02 obtained by Huckestein and Schweitzer (1994). 

The anomalous diffusive behavior can directly be seen in the spreading of wave packets 
(Chalker and Daniell, 1988; Huckestein and Schweitzer, 1994). While the variance R{2,t) = 
J (i^r |r p |0(r, t)^ of a wave packet (f){r,t) at the critical point grows linearly in time, which 
merely reflects the scale invariance of the system, the return probability p{t) = |0(O,t)p to 
the origin is enhanced and scales like 

A nonzero value of r] should be observable in measurements of the Coulomb drag in 
double layer systems (Shimshoni and Sondhi, 1994) and in hot-electron relaxation rates and 
electron-electron and electron-phonon scattering rates (Brandes et ai, 1994). 

The scaling of the two-particle spectral function S{q,uj) with q'^/uj allows to relate the 
correlation dimension D{2) of the local density to the correlation dimension D{2) of the 
spectral measure associated with the local density (Ketzmerick et ai, 1992; Huckestein and 
Schweitzer, 1994). In two dimensions, D{2) = D{2)/2. 

B. Universality of Multifractal Spectra 

In this section we review some numerical results about f{a) spectra for different models 
of the quantum Hall system and different physical quantities. We start with our prototype 



measure, the local density |'?/'(r)p. The /(«) spectrum in Fig. ^is calculated for the density 
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shown in Fig. |2^. In addition to the numerical data the parabohc approximation ( p.l9|) is 



shown with oq = 2.29 ± 0.02 (Huckestein et ai, 1992). The variation of ao for different 
eigenstates in the energy range where ^{E) ^ L is of the same order. This shows that 
the /(a) spectrum not only describes a single eigenstate but is a universal feature of the 
system studied. While these calculations were performed for a tight-binding Hamiltonian, 
the results agree with calculations by Pook and JanBen (1991) for a real space potential on a 
square of linear size 70/c who obtained ao = 2.3±0.07, Doc = 0.95±0.1, and D^oo = 3.7±0.1 
(JanBen, 1994). 

These calculations give a strong indication that the /(a) spectrum of the local density 
is universal at the critical point and that exponents, such as ao, are suited to describe the 
phase transition. It is not clear to what extend this universality holds for other physical 
observables. Calculations by Huckestein and Schweitzer (1992a) for the equilibrium current 
density |j(r)| as the local observable (see Fig. ^7\), by Schweitzer (1992) for the local mag- 
netization |m(r)| = |r X j|, and by Fastenrath, JanBen, and Pook (1992) for the Thouless 
numbers give the same spectrum within the error bars and support the notion of universality. 
On the other hand, it is easy to construct multifractals with different /(a) spectra. Given a 
normalized multifractal measure p the normalized moments p"^ generate new measures with 
a[p'"l(g) = ma{mq) - r(m) and /[p'"](a[p"]) = f{a{qm)) (JanBen, 1994). While this /(a) 
spectrum is simply related to the spectrum of p it is not identical. 

C. Multifractality and Scaling 

We now want to relate the results about the multifractal properties at the quantum Hall 
transition to the discussion about scaling in the previous sections. 

Bauer, Chang, and Skinner (1990) considered the scaling of generalized inverse partici- 
pation ratios [cf. Eq. ( p.9| )] 

p(9) = [ (iV|V^(r)|2'? (9.26) 
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as a function of system size L and distance from the critical energy AE". Wegner (1980) 
introduced critical exponents 7r(g) defined by 

P(^) oc AE^(«). (9.27) 

Using a standard finite-size scaling ansatz for P^''\ 

P^''\L,E) oc L-^^'i^Fg{L^/''AE), (9.28) 

the condition that P^^^ is finite for E ^ E^. leads to the scaling relation 

7r(g) = VT{q). (9.29) 

In particular, for the inverse participation ratio one finds 7r(2) = uD{2). Hikami (1986) cal- 
culated the exponent 7r(2) = 3.8±0.4 for the lowest Landau level. Neglecting the multifractal 
structure of the wavefunctions and assuming P^"^^ oc he concluded u = tx{2)/2 = 1.9±0.2. 
Using D{2) = 1.62 ±0.02 obtained from numerical calculations (Huckestein and Schweitzer, 
1994) and the correct relation ( |9.29| ), we get 



i/ = 2.4±0.3, (9.30) 

in agreement with the result ( p.4| ) obtained from finite-size scaling of the localization length. 

In the spirit of Shapiro (1987) one can try to study the scaling of the whole distribution 
of box-probabilities P{p, A). At the critical point this distribution becomes universal and is 
determined by the /(a) spectrum. The bulk of the distribution is log-normal with deviations 
in the tails of the distribution that account for the differences between the parabolic ap- 



proximation and the true f{a) curve in Fig. ^ In this limit the parameter characterizes 



the whole distribution. Away from the critical energy, as long as C,{E) ^ L, the distribu- 
tion will still show power-law scaling ( |9.18|) . Since the system moves towards the localized 



regime the typical box-probability ptyp "will decrease and by Eq. ( |9.21|) will increase. If 
single-parameter scaling of the distribution holds, ao has to obey a finite-size scaling law, 

ao(^) = A{L^/''AE) = ao{E,) + aL^/^^l AP| + ■ ■ ■ . (9.31) 
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The linear energy dependence is compatible with the numerical simulations presented in 
Fig. ^ (Huckestein and Schweitzer, 1992b), but these calculations do not rule out a quadratic 
energy dependence. For larger distances from the critical energy, when ^ (E) becomes smaller 
than the system size L, power-law scaling breaks down on length scales larger than ^{E). 
Eq. ( |9.18| ) no longer describes the whole distribution P{p, \) and ao loses its meaning as the 



single parameter describing the distribution. 

Employing conformal invariance in two-dimensional systems Cardy (1984) was able to 
show that the scaling amplitude = MmM^oo M / Xm of the inverse correlation length 
of the correlation function of an operator on a cylinder is related to the scaling dimension x 
of the operator in the plane by 

A"^ = 27rx. (9.32) 

For the local density one is tempted to identify Ac with the amplitude derived from the 
localization length on the cylinder and x with the exponent r]/2. This argument neglects 
the broad distribution of the local density and the corresponding broad distribution of the 
correlation functions. Instead of mapping a single correlation function from the plane onto 
the cylinder, in the presence of broad distributions one needs to map the whole distribution 
function. Since the scaling amplitude Ac was extracted from averaging the logarithm of the 
density correlation function the corresponding scaling dimension describes the scaling of the 
typical correlations x = ao — d (JanBen, 1994). Noting that, by convention, the localization 
length is defined by the absolute value of the Green function, we get 

which yields Ac = 1.1 ±0.1 in accordance with the finite-size scaling result. A similar relation 
has been derived by Ludwig (1990) for a random 2D ferromagnet. 

The discussion of multifractals can be extended to the description of non-normalized 
observables Q (Ludwig, 1990; Duplantier and Ludwig, 1991; Pook and JanBen, 1991; JanBen, 
1994). If the normalized observables pt*^! = Q/Y^f Q form a multifractal measure the 
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distribution of Q can be described by the f^^\a^^^) spectrum of pt'^l and the normahzation 
exponent X^^\ defined by 

^Q(/)ocA^"". (9.34) 

i 

Analogous to Eq. ( |9.21D the scahng of the typical value Qtyp is then given by (Janfien, 1994) 

Qtyp(0 = exp(ln(g(/))) oc (9.35) 



JanBen (1994) interpreted Qtyp(0 the value of Q for a system of size / and the average 
in Eq. (|9.35|) as an average over an ensemble of systems. If single-parameter finite-size scaling 
holds for Qtyp and Qtyp is scale invariant at the critical point then close to the critical point 

QtyAl)=Qlyv + (^l""t- (9-36) 

The box-observable Qtyp = Qtyp — Qlyp thus scales like 

Qtyp OC l'^" (9.37) 

close to the critical point. Identifying Eqs. ( |9.35| ) and ( |9.37| ) relates the critical exponent u 
of the localization length to the multifractal exponents ajj^^ and X^^\ 

l/zy = af^l+X[«l. (9.38) 

Fastenrath, JanBen, and Pook (1992) calculated Thouless numbers near the center of the 
lowest Landau level for systems of varying size. They found that the typical Thouless 
number gtyp{l) = exp[(ln (/(/))] scales according to Eq. ( |9.35| ) with a|f' = 2.25 ± 0.05 and 



X^a] = _i.75 ± 0.05. Using Eq. ( WM , v = 2.2± 0.3 in accordance with Eq. 



X. DISCUSSION AND CONCLUSIONS 

Having reviewed the experimental and theoretical results we are now in a position to 
compare both. The most important result of the comparison is that both experiments and 
numerical simulations show similar scaling behavior. The plateau transition in the integer 
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QHE can thus be understood as a continuous phase transition of noninteracting electrons 
with a single diverging length C,- If we focus on spin-split Landau levels there is a remarkable 
agreement in the value of the exponent z/, describing the divergence of the localization length 
^, between experiment [2.3 ± 0.1 (Koch et al, 1991b), 2.4 ± 0.1 (Wang et al, 1994a)] and 
numerical simulation [2.35 ± 0.03 (Huckestein, 1992)]. This suggests that the temperatures 
in these experiments are in fact low enough to observe the critical behavior and influences 
of the finite width of the Fermi-Dirac distribution are negligible. 

From the value of u and the value of l/uz = 0.41 ± 0.04 measured in the dynamic scal- 
ing experiment of Engel et al. (1993) we get the dynamic critical exponent z = 1.0 ± 0.1. 
Dynamical scaling can also be used to explain the temperature dependence of the transi- 
tion. In the approach outlined in Sec. P^I A| the temperature sets the effective system size 
Lcs oc T~P/^ that has to be compared to the localization length ^ oc \AE\~'^. In this 
interpretation the measured temperature exponent k, = p/2v. In the dynamic scaling inter- 
pretation the temperature enters as the frequency range uj = ksT /h over which excitations 
are present in the system. Scaling as a function of uot with t oc C,^ leads to the identification 
K, = l/i^z, in agreement with the results of Wei et al. (1988; 1990) from experiments on 
InGaAs/InP heterostructures and on AlGaAs/GaAs heterostructures (1992). Koch et al. 
(1991a; 1991b), using AlGaAs/GaAs heterostructures, however, find values of k, which vary 
between samples. Wei et al. (1992) argue that there is a characteristic temperature below 
which universal scaling behavior is observed and that this temperature is much lower in 
AlGaAs/GaAs than in InGaAs/InP. The agreement between the dynamic scaling experi- 
ment and the temperature scaling below the characteristic temperature suggests that these 
experiments are indeed probing the asymptotic low-temperature regime described by the 
dynamic exponent z. Which interpretation describes the experiments correctly presumably 
depends on the relative size of Les/^ and ksTr/h. 

We then need to understand why other experiments see power-law scaling with tem- 
perature but with nonuniversal exponents. This difference could be due to the different 
ranges of the disorder potentials in the different materials. In InGaAs/InP the electron gas 
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is formed in the InGaAs, which is an alloy with a potential presumably varying on an atomic 
length scale. By contrast, in AlGaAs/GaAs heterostructures the electron gas is situated in 
the GaAs, which forms a perfect crystal. In conventional heterostructures, doped in the 
AlGaAs to increase the mobility, the major scattering mechanism is remote ionized donor 
scattering. The samples used in the experiments by Koch et al. were additionally doped 
in the plane of the electron gas. Since this introduces charged impurities into the quasi- 
two-dimensional gas it is not obvious that the range of the scatterers is actually sufficiently 
different in both materials. An increase in the correlation length of the disorder potential 
leads to an increase in the irrelevant length scale ^irr and hence to an increase in finite-size 
effects. 

The temperature dependence of different phase-breaking mechanism could also account 
for different characteristic temperatures. At sufficiently low temperatures phase-breaking 
occurs presumably due to a single process. Coulomb interaction between the electrons, as- 
sociated with the dynamical exponent z. At higher temperatures other processes, like the 
electron-phonon interaction, become important. Since these processes will, in general, have 
different temperature dependences the observed nonuniversality of k might be a crossover 
phenomenon. The characteristic temperature of the crossover to the asymptotic, universal 
low-temperature regime is nonuniversal and might show a strong sample dependence. 

In addition to the relevant length scale ^ the numerical simulations show the influence 
of an irrelevant length scale For higher Landau levels and short correlation length of 
the disorder the irrelevant length can become orders of magnitude larger than the magnetic 
length and hence comparable to the phase coherence length in experiments. On length scales 
smaller than ^irr the scaling is described by a two-parameter flow according to Eq. ( ^.13|) 
while the asymptotic scaling behavior for large length scales is given by a single-parameter 
scaling relation. In a a^x vs. a^y flow diagram we expect the renormalization group flow 
to produce a set of different curves in the former regime and a single scaling curve in the 
latter. The differences in the flow diagrams of Wei et al. (1987), of Yamane et al. (1989), 
of Kravchenko et al. (1990), and of McEuen et al. (1990) might have their origin in the 
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different size of the crossover length scale. 

The value z = 1 extracted from the experiments does not agree with the result 2; = 2 for 
noninteracting electrons. This means that another mechanism is necessary to describe the 
dynamics of the system. It has been argued by Polyakov and Shklovskii (1993) and by Lee, 
Wang, and Kivelson (1993) that the Coulomb interaction between the electrons gives rise to 
z — 1. This can be motivated by the naive scaling argument that the relevant energy scale 
is the Coulomb energy on the length scale ^, E ex e^/^, which immediately provides z — 1. 
If Coulomb interactions have to be taken into account to understand the dynamic scaling 
behavior, this raises the question whether these interactions are relevant for the static scaling 
exponent too. The agreement between experiments and simulations for noninteracting 
electrons suggests that u is not changed by Coulomb interactions. However, no calculation 
of v for interacting electrons under quantum Hall conditions has been performed so far, 
although Lee, Wang, and Kivelson (1993) claim that the edge states of interacting quantum 
Hall droplets behave hke noninteracting particles and that the results of the network model 
apply in presence of interactions as well. 

For spin-degenerate Landau levels the measured values of u do not agree with the nu- 
merical value u — 2.35 (Koch et al., 1991b). Also the measured value of k is much smaller in 
the spin-degenerate case (Wei et ai, 1988; Koch et ai, 1991a; Engel et ai, 1993). The data 
analysis in all experimental investigations assumed a single critical energy in the center of 
the spin-degenerate Landau level. Polyakov and Shklovskii (1993), Lee and Chalker (1994), 
and Wang, Lee, and Wen (1994b) argued recently for the existence of two separate critical 
energies even in the completely spin-degenerate case. Lee and Chalker (1994), and Wang, 
Lee, and Wen (1994b) supported their arguments with numerical results that were compat- 
ible with two critical energies and u = 2.35. When fitting their data with a single critical 
energy, they observed a divergence of ^ with an exponent ~ 5.8 in some intermediate range 
of energies. It has to be seen whether experimental data can also be fitted with two critical 
energies and u = 2.35. 

As far as the value of the conductivity tensor at the transition is concerned there is 
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not much agreement between experiment and theory. While theoretical arguments and 
calculations seem to favor a^x = O.be^/h and a^y = (n + l/2)e^//i (at least in the lowest 
Landau level), the experimental situation is completely unclear. 

So far, there exists no direct experimental evidence for the multifractal structure of 
eigenfunctions at the transition. Since the multifractality leads to an anomalous behavior 
of the density-response function signatures of multifractality might be observable in the 
temperature dependence of transport processes (Chalker and Daniell, 1988; Brandes et al, 
1994). 

In summary, we have shown that numerical simulations are a suitable instrument to 
calculate the critical behavior of noninteracting electrons in the quantum Hall system. The 
observed behavior agrees with the description of the plateau transitions as continuous phase 
transitions with a single diverging length scale. The numerically determined critical ex- 
ponent ly — 2.35 ± 0.03 agrees well with experimental results. In addition to the relevant 
scaling field driving the transition an irrelevant scaling field is observed in the numerical 
simulations. The associated crossover length can become macroscopic, in which case the 
irrelevant scahng field has to be included in the scaling fiow. 
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APPENDIX: RANDOM-LANDAU-MATRIX-MODEL 



The generalization of the result for the correlation function Eq. ( 3.24|) to arbitrary Landau 
level index n was given by Mieck (1993), 

{nk,\V\nh){nh\V\nh) = ^ ^ exp ^-^{k, - k.ffj') exp [-\{k, - k.fllj^ 



fc=0 m=0 /=0 



,-1/2 / (fci - k,)Hl \ ( ih - k,)Hl \ ( l {k,-k,fll \ .... 

^^n-k \ r, ]^n-m\ ] ^k+m-l \ 2 /^^(/^a _ p I' ^^^l 



It is instructive to study the dependence of the correlation function on the correlation 
length of the disorder potential and the Landau level index. As the correlation length and 
hence (3 increases, the off-diagonal matrix elements become smaller but the range of the 
correlations along the diagonals increases. In higher Landau levels additional polynomials 
enter the correlation function due to the presence of the Hermite polynomials in the 
harmonic oscillator functions. While in the lowest Landau level the correlations (|3.24|) are 
always positive, in higher Landau levels the correlations start to oscillate for short correlation 
lengths and small differences ki — k^. In the limit of large cr/lc, the correlations become 
independent of the Landau level and take on the universal form given by the expression 
for the lowest Landau level, Eq. ( |3.24| ). This limit corresponds to the semi classical limit 



discussed in Section |1II D| that is independent of the Landau level index. 

In order to use the Ansatz (|3.28| ) for the matrix elements {niki\V\n2k2) the weight 
function h{x, k) appearing in it has to be expressed through the correlation function g{x, k) of 
the disorder potential. Using the translational invariance of the disorder-potential correlation 
function and the fact that the disorder potential is real, the weight function h{x, k) is seen 
to be a real symmetric function, 

h{x, k) = h{—x, k) = h{x, —k) = h{—x, —k). (A2) 

With Eq. ( |3.21| ), h{x^ k) is related to the correlation function g{x^ k) by 
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g{xi — X2, ^) — J d-^' h{xi — x', k)h{x2 — x' , k), (A3) 

which can easily be solved by Fourier transformation to give 

~g{q,k) = [hiq,k)]\ (A4) 

where g{q, k) and h{q, k) are the Fourier transforms with respect to x of g{x, k) and h{x, k), 
respectively. For the simple case of Gaussian correlations we find that h{x, k) is given by 

h{x, k) = ^^e-'=^'^^/^e-^^/-^ (A5) 
w 27rLj,cr 

and Eq. becomes 

{niki\V\n2k2) = , exp 2^ ^'ni]n2''n'^ 

'27rLy(r V 4 / V 



which gives Eq. (|3.30|) for ni = ^2 = 0. 

For numerical calculations, the integration in Eq. ( |3.30| ) needs to be discretized. Huck- 
estein and Kramer (1989) chose a step width of AX/2/c, where AX = Akl"^ = 27rZ^/Ly is 
the difference in the center coordinate of neighboring Landau wavef unctions. The discrete 
version of Eq. ( p.30| ) then reads 

(0fci|V|0fc2) = — ^exp' 



27Tlf,LyTl 



1/2 ^""i" I 4 



Y,U2k^+^M-kl exp l^-JJ^^J ' i^^) 



with 



The limits on the summation over i can be chosen such that the influence of the neglected 
terms is less than the statistical fluctuation due to disorder. 
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It remains to specify the distribution function of the random variables Ux,k- To corre- 
spond to a Gaussian distribution in real space they should also be taken from a Gaussian 
distribution. Numerically, it is faster to generate random numbers with a uniform distri- 
bution. While the distribution function of the matrix elements seems to strongly influence 
finite-size effects for the three-dimensional Anderson model (Kramer et al, 1990; MacKin- 
non, 1994), in the present case no noticeable differences could be observed between Gaussian 
and uniform distributions of the same variance. 

A similar approach to that presented above was developed by Weidenmiiller and Mieck 
(Weidenmiiller, 1987; Mieck, 1990; Mieck and Weidenmiiller, 1991; Mieck, 1993). Central 
to their idea is the observation that operators projected onto a single Landau level can 
be expanded in matrix polynomials that are related to the magnetic translation operators 
(Weidenmiiller, 1987). For the correlation function of the disorder potential this expansion 
proceeds as follows. Eq. (|3.19|) can be rewritten in terms of the Fourier transform ^(p) of 



9{r - r'),0 

/(Pp 

x(niA;i| exp[ip ■ ri] 1722/^2) ("n-s^s I exp[-ip ■ r2]|n4A;4). (A9) 

The matrix elements of the plane waves can be expressed using the magnetic translation op- 
erators. These are generated by the operator V that complements the canonical momentum 
n = p — eA. In the gauge A = -B(— (1 — a)yex + axey) these two operators are given by 

^^fp. + eBil-a)y\ p. - eBay \ 

\ Py — eBax J \py + eB{l — a)x J 

The commutation relations for these operators are those of the two-dimensional harmonic 

oscillator where the canonical momentum 11 is diagonal in the degenerate, internal states k 

of each Landau level while V is diagonal between the Landau bands n, 

[n„ n,] = leBh, [p,, Vy] = -teBh, [n, v] = o. (Aii) 



^For simplicity, we discuss the limit Ly ^ 00 in the following. 
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With these operators ip ■ r can be expressed as 



2p-r = 4^(e, xp)(P-n). (A12) 
en 



Using the commutation relations Eq. ( |A11| ) this can be exponentiated to give 



exp[zp ■ r] = exp[{i/eB){ez x p){V — 11)] 

= exp[(i/e5)(e^ x p)V] exp[-{i/eB){e, x p)n]. (A13) 

The matrix elements of exp[zp ■ r] are then given by the matrix elements Ukj^^k-^ip) = 
{niki\exp[i/eB{e^ x p)V]\n2k2) and PVni.nalp) = {niki\exp[-i/eB{e^ x p)n]|n2A;2), 

(nifcil exp[ip ■ r]\n2k2) = Uk^^kiipWn^^niip)- (A14) 

Note that due to the commutation relations, t^fci,fc2(p) does not depend on the Landau- level 
index, and W^ni,n2(p) does not depend on the internal quantum numbers k. Using Eq. ( [A13D , 



the correlation function ( |A9| ) can now be expressed in terms of the matrix polynomials 

7^^(|P|) 

ZTT 

^Uk,,k2{p)Wn,,n2{p)Uk3Mi-PWn3,n4i-p)- (Al5) 

The particular form of the matrix polynomials depends on the choice of the gauge. For the 
Landau gauge ( |3.13| ) they are given by (Mieck, 1993) 

Uk,k'{p) = exp[tp,{k + k')/2]5ip., ~k + k') (A16) 



u = mm 



m{n,n), tan (p = py/px- (A18) 



The delta function in the matrix polynomial Uk^k'ip) reproduces the delta function that we 
already observed in the correlation function ( |3.21| ). 

The matrix-polynomial representation ( |A14| ) of the plane-wave operator can also be used 
to generate the matrix elements {niki\V\n2k2)- For this, we express the real space potential 
V^(r) by its Fourier transform V^(k), 



V(v) = I — e-*'^-^F(k). (A19) 

271 



Using Eq. ( |3.18|) the second moment of the Fourier transform is given by 



V{k^)V{k2) = 27r^(|ki|)5(ki + ks). (A20) 

Thus the Fourier transform of the random potential is 5-correlated and can be expressed by 
Gaussian random variables ?)(k) of unit variance, 

v{^) = {2ny/'mk\)]'^"m- (A2i) 

The potential V^(r) is then determined by ^(k), 

V{r) = J mm'^' e-^-^(k), (A22) 

where the second moment of v{k.) is given byQ 



{}*(ki)C(k2) =5(ki -ks) (A23) 

Using Eq. ( [A14| ) the matrix elements of the disorder potential can now be expressed in 
terms of the random variables ?)(k) and the matrix polynomials, 

{n,h\V\n,k2) = I mp\)f' f/n.n.(-p)l^.,..(-p)^5(p). (A24) 



^Note that since V{r) is real i;*(k) = {;(— k). 
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FIGURES 

FIG. 1. The temperature dependence of the longitudinal and Hall conductivity [from (Kawaji 
and Wakabayashi, 1987)]. 

FIG. 2. Longitudinal conductivity Oxx vs. Hall conductivity Oxy in a Si-MOSFET for temper- 
atures between 0.35 K and 1.5 K [from (Yamane et al., 1989)]. 

FIG. 3. Longitudinal conductivity axx vs. Hall conductivity axy in an InGaAs-InP heterostruc- 
ture for temperatures between 50 mK and 770 mK [from (Wei et al, 1987)]. 

FIG. 4. Hall resistivity pxy, longitudinal resistivity pxx, and the derivative dpxy/dB of an 
InGaAs-InP heterostructure [from (Wei et al., 1988)]. 

FIG. 5. The temperature dependence of the maxima of dpxy / dB and the width ^13 of the pxx 
peaks from Fig. ^ The slope of the straight lines gives {dpxy/dB)^^^ oc T^'^ and IS.B oc T'^ with 
K = 0.42 ± 0.04 [from (Wei et al, 1988)]. 

FIG. 6. Temperature dependence of the peak width /S.B for samples of different size [from 
(Koch et al, 1991a)]. 

FIG. 7. Real part of the dynamical conductivity Re((Txx) at different frequency and tempera- 
tures [from (Engel et al, 1993)]. 

FIG. 8. Regions where V{r) < E are shaded. The eigenfunctions are localized within a distance 
Ic of the edge of the regions. The arrow indicates the direction of guiding center motion [after 
(Trugman, 1983)]. 

FIG. 9. Equipotential lines of a smooth random potential. The thick lines correspond to the 
links in the network model. The circles enclose the saddle points [after (Chalker and Coddington, 
1988)]. 
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FIG. 10. The network model of Chalker and Coddington. Each node represents a saddle point 
and each link an equipotential line of the random potential [after (Chalker and Coddington, 1988)]. 



FIG. 11. The renormalization-group flow near a saddle-point fixed point xi = X2 = 0. 
FIG. 12. j3 function according to Abrahams et al. (1979). The slope of the dashed line is l/u 



[cf. Eq. iWM 



FIG. 13. Renormalization flow diagram of axx vs. axy [after (Khmel'nitskii, 1983)]. 

FIG. 14. Normalized localization length \m/M in the lowest Landau level for 5-correlated 
potential as a function of system width M (in units of \/^lc) for energies 0.01 (*), 0.05 (v), 0.07 
(a), 0.1 (□), 0.18 (o), 0.30 (★), 0.5 (o), and 1.0 (•) (in units of F) [after (Huckestein and Kramer, 
1990)]. 



FIG. 15. The scaling function \m{E)/M = A{M/^{E)) [from (Huckestein, 1992)]. 

FIG. 16. The localization length ^{E) as a function of energy normalized by the band width 
Fcr for n = 0, a = (v), n = 0, a = Ic (a), and n = 1, a = Ic (□) [after (Huckestein, 1992)]. 

FIG. 17. Normalized localization length X^j/M in the lowest Landau level for correlated po- 
tential (cr = Ic) as a function of system width M (in units of \/^lc) for energies (*), 0.03 (x), 
0.05 (v), 0.07 (a), 0.1 (□), 0.14 (*), 0.18 (o), 0.3 (o), and 0.5 (•) (in units of F) [after (Huckestein 
et al, 1992)]. 

FIG. 18. Normalized localization length Xm/M in the second (re = 1) Landau level for corre- 
lated potential (ex = Ic) as a function of system width M (in units of \/27r/c) for energies 0.01 (*), 
0.03 (v), 0.05 (a), 0.07 (□), 0.1 (o), 0.18 (★), 0.3 (o), and 0.5 (•) (in units of F) [from (Huckestein, 
1992)]. 
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FIG. 19. Normalized localization length \m/M in the second (n = 1) Landau level for uncor- 
related potential {a = 0) as a function of system width M (in units of \/27r/c) for energies 0.01 (*), 
0.1 (v), 0.18 (a), 0.3 (□), 0.5 (o), 0.65 (★), 0.8 (o), and 1.0 (•) (in units of T) [from (Huckestein, 
1992)]. 

FIG. 20. The normalized localization length \m/M at the band center of the second Landau 
level as function of /^^ = (ci^ + for system width M = 16 (o), M = 32 (•), M = 64 (o), and 

M = 128 (*) [after (Huckestein, 1992)]. 

FIG. 21. Scaling towards the fixed point value Ac = 1.14 as a function of M/^irr for the data 
shown in Fig. |2^ [from (Huckestein, 1994)]. 

FIG. 22. The irrelevant length scale (in units of \/^lc) as a function of [from (Huckestein, 
1994)]. 

FIG. 23. Deviations of the scaling function A from its asymptotic value Ac [from (Huckestein, 
1994)]. 

FIG. 24. Phase diagram of the lowest Landau level in the presence of a periodic potential with 
a = 3/5 unit cells per flux quantum. Indicated are the positions of the critical energies and the 
Hall conductivity in the different phases. The bars at zero T/Eq show the subbands in the absence 
of disorder. 

FIG. 25. The local density |V(r)|^ of an eigenstate at the center of the lowest Landau level of 
a lattice system with 150 x 150 sites [from (Huckestein and Schweitzer, 1992a)]. 

FIG. 26. The /(a) spectrum of the local density of the eigenstate shown in Fig. |2^. Crosses 
show the errors in f[a{q)) and a{q) for q = —2, -1, 1, and 2 (from right to left) [from (Huckestein 
et al, 1992)]. 
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FIG. 27. The f{a) spectrum of the local equilibrium current density of the eigenstate shown 
in Fig. |2^. Crosses show the errors in f{a{q)) and a{q) for q = —2, -1, 1, and 2 (from right to left) 
[from (Huckestein and Schweitzer, 1992a)]. 

FIG. 28. The energy dependence of ao{E) near the center of the lowest Landau band. The 
straight line is a guide to the eye [from (Huckestein and Schweitzer, 1992b)]. 
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